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COMPUTER  SIMULATION  MODEL 
of  the 
IPSWICH  RIVER  BASIN 

by 

Alan  W.  Burns  &  Ivan  C.  James  II 

ABSTRACT 


Consideration  of  all  proposals  made  by  towns  and  regional  planners 
for  development  of  water  resources  in  the  Ipswich  River  basin  is  impracti- 
cal using  traditional  hydrologic  techniques.   A  computer  simulation  model 
has  been  written  to  assess  the  different  development  proposals.  The  model 
simulates  the  hydrology  and  operation  of  water-supply  systems  in  the  basin 
at  given  levels  of  demand  and  particular  sequences  of  streamflows.  Hydro- 
logic  simulation  distributes  incremental  gaged  runoff  to  each  tributary 
drainage  area  and  computes  the  effects  of  groundwater  withdrawals  on 
streamflow.   Operational  simulation  of  the  water-supply  systems  takes 
water  from  supply  sources,  diverts  streamflow  to  storage  reservoirs,  and 
transfers  water  between  reservoirs.   Established  operating  rules  are 
observed;  for  example,  no  diversions  are  allowed  below  a  minimum  discharge, 
and  diversions  are  allowed  only  in  winter  and  spring.  Model  results  Indi- 
cate the  relative  abilities  of  different  developmental  plans  and  provides 
a  rational  means  of  decision  making  to  be  used  with  costs,  legislative 
feasibility,  and  environmental  impact. 


INTRODUCTION 


Communities  in  and  near  the  Ipswich  River  basin  are  facing  water- 
supply  problems  caused  by  the  rapid  population  growth  since  1950.   The 
expected  expansion  of  the  Boston  metropolitan  area  will  accentuate  these 
problems.   Even  the  adequacy  of  the  water-supply  systems  of  towns  whose 
population  has  remained  level  has  been  reduced  by  increased  per  capita 
consumption  and  increased  utilization  of  the  same  basic  water  resource, 
the  Ipswich  River,  by  other  towns  upstream.   Plans  for  development  of 
additional  supply  sources  are  being  considered  in  several  towns  and  at 
the  regional  and  state  level.   Some  of  these  plans  may  conflict  with  other 
alternatives  for  integrated  management  of  the  basin's  resources.   The 
well-established  need  for  additional  sources,  along  with  the  substantial 
lead  time  for  authorization,  financing,  and  construction  have  led  to 
several  regional  proposals. 
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A  necessary  requirement  for  effective  development  of  water-supply 
sources  in  this  basin  is  a  method  for  evaluating  the  various  alternatives 
that  have  been  proposed.   To  meet  this  requirement,  a  computer  simulation 
model  was  written  to  assess  alternaitve  plans  for  development  of  water 
resources  within  the  basin. 

The  model  simulates  the  operations  of  each  town's  supply  system  in 
meeting  its  demand  with  the  use  of  surface  water  and/or  groundwater.   The 
model  evaluates  sets  of  supply  systems  at  desired  levels  of  demand. 
Results  obtained  from  the  model  indicate  success  or  failure  for  a  given 
set  of  supply  systems  and  level  of  demand  -  a  failure  being  a  water  deficit 
for  one  or  more  towns.   Thus,  the  model  can  be  used  by  water  resource 
manages  to  compare  and  judge  the  effectiveness  of  alternatives.   Although 
the  model  does  not  specify  the  best  alternative,  it  provides  valuable 
information  for  the  manager,  who  must  then  consider  other  factors,  such 
as  engineering  efficiency,  costs,  legislative  actions,  financing,  inter- 
town  cooperation,  pollution,  and  environmental  impacts. 

The  model  has  been  designed  for  the  Ipswich  River  basin,  but  lias  been 
kept  as  general  as  possible  for  use  in  other  basins.   Three  physical 
characteristics  of  the  Ipswich  River  basin  that  must  be  present  in  another 
basin  to  use  the  model  are: 


1)  one  main  artery  -  because  no  major  tributaries  enter  the 
Ipswich  River,  the  model  was  developed  without  branching 
routines.   However,  the  modification  needed  to  model  a 
basin  with  more  than  one  main  branch  would  not  be  diffi- 
cult to  program. 

2)  localized  aquifer  -  because  most  of  the  groundwater  is  in 
small,  isolated  aquifers  with  close  hydraulic  connection 
to  the  stream,  the  model  was  not  written  to  simulate  the 
entire  aquifer  but  to  simulate  the  response  of  streamflow 
to  pumping  wells.  More  complex  models  are  required  in 
basins  with  extensive  aquifers. 

3)  off-channel  surface-water  storage  -  because  there  are  no 
adequate  sites  for  on-channel  reservoirs,  the  model  was 
developed  to  use  only  pumped  storage  reservoirs.   However, 
the  modification  needed  to  model  a  basin  with  mainstream 
reservoirs  would  not  be  difficult  to  program. 


METHODOLOGY  OF  THE  MODEL 


Parameters  describing  the  physical  characteristics  of  the  supply- 
systems,  the  disposal  systems,  and  the  river  basin  are  data  required  for 
the  model.   Physical  objects,  such  as  tributaries,  reservoirs,  points  of 
diversion,  wells,  gaging  station,  and  towns,  are  coded  in  a  spatial 
concept,  which  allows  the  computer  to  interpret  geographical  locations. 
The  river  basin  is  divided  into  several  smaller  subbasins,  each  identified 
by  a  point  on  the  river,  which  is  referred  to  as  a  node  (fig.  1).   These 
subdivisions  are  made  so  that  every  physical  object  of  the  same  type  must 
lie  within  different  subbasins.   For  example,  two  reservoirs  must  lie  within 
two  separate  subbasins.   However,  a  subbasin  may  contain  all  the  different 
types  of  physical  objects.   For  example,  a  town,  a  reservoir,  a  well,  and 
a  tributary  may  lie  within  the  same  subbasin  (fig.  2).   The  exception  to 
this  is  that  any  number  of  wells,  whose  effects  are  on  the  same  tributary, 
may  be  within  the  same  subbasin.   Wells,  reservoirs,  or  towns  that  are 
outside  the  drainage  basin  of  the  river  but  that  still  affect  the  basin  in 
some  way  are  coded  within  subbasins  defined  to  have  zero  drainage  area 
(relative  to  the  river  being  modeled)  (fig.  2).  After  all  subbasins  have 
been  defined,  their  respective  nodes  are  numbered  in  downstream  order. 


^^  Tributary  Node 
O    Reservoir  Node 
^f  Diversion  Node 
#  Well  Node 
(•)  Town  Node 
A  Gaging  Node 


Figure   1 


SUBBASINS  IDENTIFIED   BY   DIFFERENT  TYPES  OF   NODES 


MM  Sub-basin  +  Tributary  &  Reservoir 

WM  Sub-basin  +  Out-of-basin  Reservoir  &  Town 

11  Sub-basin  +  Town 


o 

0 


Sub-basin  +  Out-of-basin  Weli  &  Town 
Sub-basin  +  Well 


Figure  2 


MULTIPLE  USE  OF  NODES 


The  model  simulates  the  continuous  interactions  of  the  real  world 
by  cycling  (fig.  3)  through  discrete  "passes"  (a  pass  being  a  specific 
function  executed  on  a  subset  of  nodes  in  downstream  order)  on  a  monthly 
time  interval.  Monthly  flows  for  each  of  the  two  gaging  stations  within 
this  basin  are  the  hydrologic  input  for  the  model. 

The  first  pass  distributes  the  incremental  runoff  of  the  gaging 
stations  to  each  subbasin. 


Rg  =  Qg/DAg 


where  R  =  incremental  runoff 

Q  =  incremental  discharge 
DA=  incrimental  drainage  area 
g  =  subscript  indicating  gaging  station  index 


Qi     =  Rj  x  DAi     where    i  =  subscript  indicating  nodal  index 

j  =  subscript  indicating  respective  gaging  station  index 
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Figure  3 


FLOWCHART  OF  THE  MODEL 


When  there  is  a  reservoir  within  a  subbasin,  the  monthly  inflow  is 
computed  from  the  drainage  area  of  that  reservoir  and  the  average  of  the 
runoff  for  all  the  discharge  records. 


I^  =  R  x  DAk  x  C  where  I  =  reservoir  increase 

ft  =  average  of  the  runoff  for  all  the  discharge  records  used 
C  =  constant  to  convert  rate  to  volume  per  month 
k=  subscript  indicating  reservoir  index 


If  water  is  diverted  from  another  river  to  augment  flow,  it  is  added 
in  at  the  node  corresponding  to  the  real  geographic  location. 

The  second  pass  accounts  for  streamflow  depletion  resulting  from 
groundwater  withdrawals.   Streamflow  depletion  computed  by  a  method  des- 
cribed by  Jenkins  (1970)  is  based  on  well  pumpage,  distance  from  the  well 
to  the  stream,  transmissivity ,  and  specific  yield.   The  cumulative  effect 
on  stream  flow  is  obtained  by  summing  the  effects  of  each  month's  pumpage. 
Whether  this  depletion  of  streamflow  is  interpreted  as  interception  of 
groundwater  that  would  have  become  streamflow  or  as  infiltration  of  exist- 
ing streamflow  is  not  important  to  the  model. 

Characteristically,  the  aquifers  are  recharged  every  spring  and  reach 
their  maximum  yearly  storage  during  spring  floods  (Sammel  and  others  1966)  . 
To  simulate  the  spring  recharge,  the  accumulated  effects  on  streamflow 
for  each  well  are  reduced  to  zero,  and  an  equal  amount  of  surface  runoff 
is  prevented  from  becoming  streamflow. 

The  third  pass  attempts  to  satisfy  the  demand  of  each  town  through 
use  of  a  set  of  normal  supply  sources,  which  may  include  wells,  streamflow 
diversions,  and  surface  reservoir  withdrawals.   Each  town  draws  water  from 
its  normal  supply  sources  in  accordance  with  its  list  of  priorities  until 
the  demand  is  satisfied  or  until  the  sources  are  depleted  or  tiieir  capaci- 
ties exceeded.  Any  unsatisfied  demand  is  carried  over  to  the  next  pass  for 
further  operation.   As  water  is  used  by  each  town,  a  fraction  of  it  returns 
to  the  streamflow  according  to  the  value  of  its  sewage  function  for  that 
month.   (See  System  Parameter  Definition.) 

The  fourth  pass  attempts  to  satisfy  any  demand  not  met  by  the  normal 
supply  sources  from  a  set  of  extended  supply  sources.   These  sources  may 
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include  wells,  surface  reservoir  withdrawals,  and  purchase  of  water  from 
a  source  having  no  effect  on  the  basin's  hydrology.   Any  town  that  was 
unable  to  meet  its  demand  through  its  normal  supply  sources  attempts  to 
fulfill  the  remaining  demand  from  its  extended  sources  until  the  demand  is 
met  or  the  sources  are  depleted  or  their  capacities  exceeded.  Any  part  of 
the  demand  still  not  satisfied  is  recorded  as  a  failure  of  the  system  for 
that  month. 

The  fifth  pass  diverts  streamflow  to  pumped  storage  reservoirs.   Three 
factors  limit  the  amount  of  water  diverted  for  each  reservoir: 

1)  The  amount  of  water  needed  to  fill  the  reservoir. 

2)  The  constricting  factor  of  pump  and  pipe  size. 

3)  The  amount  of  water  above  the  legally  defined  minimum 
discharge  that  must  flow  past  the  diversion  point. 

The  sixth  pass  diverts  water  from  one  reservoir  to  another.   This  is 
used  for  several  systems  where  the  stream  diversion  is  to  a  reservoir  that 
is  not  directly  drawn  on  by  the  treatment  plant.   Three  factors  limit  the 
amount  of  water  diverted  from  one  reservoir  to  another: 

1)  The  amount  of  water  needed  to  fill  the  receiving  reservoir. 

2)  The  constricting  factor  of  pump  and  pipe  size. 

3)  The  amount  of  water  in  the  distributing  reservoir. 

Upon  completion  of  these  six  passes,  the  monthly  time  interval  is 
incremented  to  the  next  month,  a  new  set  of  mean  monthly  flows  are  read 
into  the  model,  and  the  entire  set  of  operating  passes  is  repeated.   This 
cycling  is  continued  for  the  prescribed  number  of  months  in  the  simulation 
period. 

SYSTEM  PARAMETER  DEFINITION 


Extensive  data  collection  is  necessary  to  provide  physical  and  opera- 
tional parameters  for  the  simulation  program.   The  model  structure  and  the 
level  of  detail  required  in  the  data  are  designed  to  make  the  results  fit 
the  scale  of  decisions  faced.   If  every  detail  were  used,  the  entire  pro- 
ject would  be  merely  a  data  collecting  project.   The  model  would  be  unique 
and  could  not  be  used  in  another  basin,  and  its  usefulness  in  decision 
making  would  be  lost  in  the  complexity  of  details.   On  the  other  hand,  the 
data  must  be  accurate  and  specific  enough  to  describe  the  pertinent  opera- 
ting characteristics  of  the  water-supply  systems.   Refinement  of  the  data 
and  parameters  are  carried  only  to  the  extent  that  further  refinements  do 
not  improve  the  results. 
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Some  of  the  parameters  required  for  modeling  are  measurable  physical 
characteristics.  These  include  the  drainage  area  of  each  subbasin  and 
gaging  station.  The  usable  storage  capacity  and  drainage  area  of  reser- 
voirs are  recorded.   Special  codes  in  the  data  input  indicate  whether  a 
reservoir  is  outside  the  basin's  drainage  boundary  and  whether  a  reservoir 
receives  diverted  river  water  or  water  transferred  from  another  reservoir. 
The  maximum  capacity  of  either  pump  sizes  or  pipe  sizes  between  river  and 
towns,  river  and  reservoirs,  two  reservoirs,  and  reservoirs  and  towns  is 
necessary.   Wells  must  be  located  within  their  subbasin  and  the  down- 
stream, noncontributing  drainage  area  computed.  The  pumping  capacity  of 
each  well  is  recorded,  and  a  special  code  indicates  which  town  uses  the 
well.   The  effects  on  streamflow  caused  by  well  pumpage  are  computed  for 
each  well  from  the  transmissivity  and  specific  yield  of  the  aquifer  and 
distance  to  the  stream.   A  special  code  indicates  the  subbasin  for  each 
town  and  each  model  to  be  used  as  an  observation  node. 

Several  other  parameters  are  not  directly  measurable  and  require 
engineering  analysis  and  hydrologic  judgment.   These  parameters,  which  are 
held  constant  and  used  statically  in  discrete  steps  in  the  model  although 
in  reality  they  are  dynamic,  may  cause  the  deviations  between  simulated 
results  and  actual  occurrences. 

A  very  important  parameter  in  the  model  is  the  demand  for  water. 
Demand  is  not  constant  throughout  the  year,  so  an  average  demand  is 
distributed  by  a  time-demand  relationship  (fig.  4).   The  curve  computed 
for  each  town  is  based  on  three  representative  years,  in  which  demand  was 
not  affected  by  water-use  restrictions.   The  model  does  not  consider  maxi- 
mum daily  demands  of  each  town  system  but  considers  these  demands  a  prob- 
lem of  the  town's  distribution  system  rather  than  the  supply  system. 
Projected  average  daily  demands  from  1975  to  2020  in  5-year  increments  for 
each  town  are  needed  for  modeling  proposed  supply  systems  for  stage- 
development  analysis.   For  years  in  which  estimates  are  not  available, 
average  daily  demands  must  be  interpolated  or  extrapolated. 

The  sewage  return  functions  for  each  town  is  an  estimate  of  the  frac- 
tion of  the  water  pumped  to  the  town  that  returns  to  streamflow  through 
sewage  or  septic  tank  effluent.   The  percentage  is  zero  for  sewered  towns 
in  which  the  sewage  is  pumped  out  of  the  basin.   For  sewage  treatment 
plants  within  the  basin,  the  return  function  is  computed  under  the  assump- 
tion of  1)  complete  return  of  water  used  by  households ,  and  2)  loss  of 
water  used  for  outside  purposes.   For  septic  tank  effluent,  a  function  was 
developed  that  returns  most  of  the  used  water  to  streamflow  during  the 
winter  and  spring  but  reduces  the  percentage  to  a  low  level  in  summer, 
reflecting  both  outdoor  use  of  water  and  increased  evapotranspiration.   The 
septic  tank  function  and  proper  sewered  function  are  weighted  according  to 
the  percentage  of  water  users  that  are  sewered. 
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Figure  4 


Another  important  parameter  is  the  operating  rules  of  a  town  in  using 
its  supply  sources.   These  rules,  which  are  governed  by  legislative  acts, 
engineering  feasibility,  and  human  tradition,  are  variable  in  time,  but  in 
the  model  they  must  be  held  constant.  At  every  diversion  point  on  the 
river,  the  legislature  has  established  minimum  discharges,  which  must  be 
exceeded  before  pumping  is  allowed.   The  legislature  restricts  the  months 
when  diversions  are  allowed.   Within  the  model,  water  is  considered  to  be 
drawn  from  only  one  source  at  a  time,  and  thus,  all  sources  are  coded  in 
order  of  priority.  This  priority  list  is  formed  to  approximate  actual 
procedures  as  closely  as  possible. 


HYDR0L0GIC  DATA  SELECTION 


Performance  of  a  water-supply  system  is  obviously  directly  related  to 
hydrologic  conditions  in  the  basin.   Because  hydrologic  conditions  are 
dynamic  and  unpredictable,  it  is  not  reasonable  to  design  a  system  from  one 
sequence  of  hydrologic  events  and  to  expect  them  to  be  repeated  in  the 
future.   For  example,  water-supply  systems  that  respond  well  to  a  long,  mild 
drought  may  become  critically  short  during  a  short,  severe  drought.   There- 
fore, the  response  of  a  system  to  one  particular  sequence  of  hydrologic 
events  may  not  be  a  particularly  good  measure  of  the  reliability  of  that 
system. 
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Severity  and  duration  of  dry  periods  must  be  considered  in  estimating 
the  frequency  of  failure  of  a  water-supply  system.  To  achieve  any  stability 
in  this  measure  of  the  system's  performance,  many  more  combinations  of 
events  are  needed  than  are  usually  available  in  historic  records.  The  costs 
of  improving  a  supply  system  so  that  it  can  survive  a  drought  of  certain 
severity  are  similar  to  purchase  of  insurance — the  higher  the  premium  paid, 
the  greater  the  coverage  of  the  policy.   Because  these  insurance  invest- 
ments in  an  improved  system  may  be  large,  a  balance  must  be  struck  between 
accepting  some  risk  and  investing  in  the  next  unit  of  a  supply  system 
before  its  need  is  fully  realized. 

The  hydrologic  input  for  the  model  is  a  set  of  natural  mean  monthly 
flows,  as  recorded  at  the  gaging  stations,  adjusted  for  upstream  diversions. 
Selection  of  the  set  of  discharges  to  be  used  influences  the  results  of  the 
model  significantly.   The  10  years  of  historic  record  from  1960  to  1969, 
used  for  most  of  the  model  runs,  includes  some  normal  years,  a  severe 
drought,  and  the  maximum  instantaneous  peak  of  record.  The  40  years  of  historic 
record  from  1931  to  1970  was  used  to  check  the  performance  of  hypothetical 
supply  systems  on  droughts  of  shorter  duration  and  severity  than  that  of 
the  mid  1960's. 

Model  results  are  influenced  by  assumed  initial  reservoir  storages. 
The  assumed  initial  storage  in  each  reservoir  should  reflect  some  normal 
or  typical  condition  for  the  supply  systems  under  consideration.   The 
method  used  to  achieve  typical  initial  storage  conditions  is  to  simulate 
2  years  of  water-supply  operation  with  median  monthly  discharge  data. 
Resulting  final  reservoir  storages  from  this  period  are  used  as  initial 
storages  for  the  period  of  simulation  for  which  performance  statistics 
are  kept. 

Use  of  synthetic  flows  can  be  advantegeous  when  a  response  of  a 
supply  system  to  a  large  number  of  possible  sequences  of  hydrologic  input 
is  needed.   These  flows  are  generated  in  a  way  that  maintains  the  pertinent 
statistics  of  the  historical  record  (Fiering  &  Jackson  1971) .   Synthetic 
flows  can  be  used  for  a  more  detailed  investigation  of  a  small  number  of 
proposed  supply  systems  after  screening  several  proposals,  using  the 
historical  hydrology. 


MODEL  RESULTS 


The  model  can  be  used  by  the  water  resource  manager  only  if  he  can 
interpret  the  computer  output  correctly.   To  interpret  the  output,  he  must 
first  understand  how  to  read  the  tables  (see  Appendix  IV  for  examples) 
and  know  what  output  options  are  available. 


12 


I 

I 
I 
I 


The  primary  table  in  the  computer  printout  lists: 

1)  nodes  along  the  river  by  some  convenient  name 

2)  usable. storage  capacity  of  a  reservoir  within  that  subbasin 

3)  maximum  storage  level  reached  by  a  reservoir  during  the 
entire  simulated  hydrologic  period 

4)  minimum  storage  level  reached  by  a  reservoir  during  the 
entire  simulated  hydrologic  period 

5)  total  amount  of  water  diverted  from  the  river  to  a  reservoir 
within  the  subbasin 

6)  total  amount  of  water  brought  in  from  outside  the  basin  by 
a  supply  system 

7)  maximum  monthly  amount  of  water  brought  in  from  outside  the 
basin  by  a  supply  system  during  the  entire  simulated  hydro- 
logic  period 

8)  accumulated  deficits  for  a  town  during  the  entire  simulated 
hydrologic  period. 

The  only  other  listing  that  is  not  optional  to  the  model,  is  an 
automatically  printed  line  of  the  town  name,  the  amount  of  a  water  deficit 
and  the  month  it  occurred,  every  time  a  supply  system  fails  to  meet  a 
demand . 

The  optional  listings  are  triggered  by  certain  codes  in  the  input 
deck.   These  include: 

1)  A  two-part  table,  in  which  the  first  part  lists  each  town  with 
its  supply  system  and  each  source  in  order  of  its  priority  and 
its  capacity;  the  second  part  lists  every  reservoir  that  either 
receives  water  diverted  from  the  river  or  that  transfers  water 
to  another  reservoir.  The  pump-pipe  size  is  also  listed. 

2)  A  printed  line  giving  the  quantity  of  water  imported  into  the 
basin  to  augment  the  river  each  month. 

3)  A  printed  line  giving  the  quantity  of  water  diverted  by  a 
reservoir  from  the  river  each  month. 

4)  A  printed  line  giving  the  quantity  of  water  transferred  from 
one  reservoir  to  another  each  month. 


13 


I 


5)  A  table  that  lists  maximum  monthly  discharge  and  minimum  monthly 
discharge  during  the  entire  simulated  hydrologic  period  and  a 
flow-duration  table  for  any  node  desired. 

6)  A  table  that  lists  the  monthly  discharge  and/or  the  storage 
level  at  any  node  desired  for  the  enire  simulated  hydrologic 
period. 

7)  A  hydrograph  of  the  monthly  discharge  and/or  storage  level 

at  any  node  desired  for  the  entire  simulated  hydrologic  period. 

8)  A  special  graph  of  any  other  time  series  from  the  model 
operation.   However,  this  option  takes  some  understanding 

of  programming,  as  it  calls  a  dummy  subroutine,  and  the  series 
to  be  graphed  must  be  defined  in  this  subroutine. 

Although  most  of  the  listings  and  tables  are  self-explanatory,  there 
are  a  few  idiosyncrasies  of  the  model.   For  instance,  a  reservoir  may  go 
dry  as  a  town  draws  water  from  it,  but,  before  the  statistics  for  that 
month  are  updated,  that  reservoir  may  have  an  inflow  from  another  reservoir 
or  from  a  diversion  of  the  river.   In  this  case,  if  the  reservoir  is  the 
town's  only  source,  a  deficit  would  have  been  printed,  and  yet  the  reserr 
voir  would  not  be  listed  as  ever  having  gone  dry.   These  apparent  ambigu- 
ities come  from  the  necessary  simplifications  made  in  the  mathematical 
concepts  of  the  model;  the  pertinent  one  in  this  case  is  the  discrete 
monthly  approximation  of  processes  that  happen  continuously.   This  approxi- 
mation emphasizes  the  importance  of  experience  and  familiarity  in  inter- 
preting the  model  results. 

Care  should  be  taken  to  check  the  frequency  and  severity  of  deficits 
in  the  results  of  the  model.   Deficits  may  reflect  a  constraint  in  the 
supply  system,  a  severe  hydrologic  condition,  or  possibly  a  mistake  in 
coding  the  intended  supply  system  for  that  run.   Each  of  these  has  quite  a 
different  meaning,  and  the  analyst  should  check  each  set  of  outputs 
carefully  before  interpreting  that  run  or  preparing  the  next  run. 
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SUMMARY 


A  simulation  model  to  aid  in  water  resource  planning  for  the  Ipswich 
River  basin  has  been  developed  and  utilized.  The  functional  operation  of 
the  model  components  and  the  procedures  for  converting  an  alternative 
development  plan  into  input  parameters  to  the  computer  model  are  described 
to  allow  the  reader  to  become  familiar  with  the  philosophy  and  capabilities 
of  the  model.   A  proficient  operating  ability  with  this  model  can  be  gained 
only  through  a  thorough  investigation  of  the  model  structure  and  some  ex- 
perience in  its  use,  based  on  documentation  provided  in  the  appendices. 

Application  of  the  model  to  any  particular  problem  is  the  responsi- 
bility of  the  user.  Reasonableness  of  the  combinations  and  linkages  of  the 
chosen  supply-system  components  are  particularly  important.   Selection  of 
criteria,  such  as  the  hydrologic  sequence  to  design  for,  may  have  a  sub- 
stantial impact  on  the  model  results  in  some  cases. 

The  model  is  now  ready  for  use  in  assessing  many  proposals  on  a 
quantitative  basis.   Constructive  criticism  of  this  report  may  lead  to 
changes  and  additions  to  the  model,  some  of  which  are  already  under  con- 
sideration. Future  developments  in  this  study  will  have  to  consider  whether 
activities  such  as  costing,  environmental  impact  assessment,  and  benefit 
distribution  should  be  done  internally  or  externally  from  the  computer 
model. 
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APPENDIX  I 


DEFINITION  OF  PROGRAM  VARIABLES 

The  following  is  a  list  of  the  more  important  variables  within  the 
program  and  their  definitions. 


A(I) 

ADD   . 

C(I) 

CLASS (I) 
CONSUP(I) 
D(I) 
DA  (I) 


-  The  plotting  value  for  one  of  A  time  series  in  the  Ith  month. 

-  The  amount  of  water  returned  to  streamflow  from  the 
sewage  function. 

-  The  plotting  value  for  one  of  A  time  series  in  the  Ith  month. 

-  The  plotting  value  for  one  of  A  time  series  in  the  Ith  month. 

-  The  Ith  discharge  used  as  an  interval  for  the  frequency  table 

-  The  average  daily  demand  for  the  Ith  town. 

-  The  plotting  value  for  one  of  A  time  series  in  the  1th  month. 

-  The  drainage  area  of  the  Ith  discharge  record  station. 


DELTA      -  The  width  of  an  interval  during  plotting. 
DEMAND (I, J)-  Two  definitions: 

1)  when  read  in,  the  percentage  of  the  average  daily  demand 
for  the  Ith  town  in  the  Jth  month. 

2)  when  multiplied  by  the  respective  average  daily  demand 
(CONSUP),  the  demand  for  the  Ith  town  in  the  Jth  month. 


DIS(I) 

DIVERT 
EFF(I,J) 
FORMT 
FREQ(I,J) 

FRMAT 

FR11ATT 


-  The  mean  monthly  discharge  of  the  1th  discharge  record  station. 

-  The  amount  of  water  diverted  from  the  river. 

-  The  accumulated  well  effects  for  the  Ith  well  in  the  Jth  month. 

-  The  values,  used  for  setting  up  variable  formats  FRMAT  and  FRMATT. 

-  The  percentage  of  time  the  discharge  of  the  Ith  node  is 
within  the  Jth  interval. 

-  A  variable  format  to  print  the  node  name  for  a  flow  trace  or 
storage  level  trace. 

-  A  variable  format  to  print  the  discharge  or  storage  level 
for  a  flow  trace  or  a  storage  level  trace. 
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GWMAX(1,J)-  The  well  pumpage  capacity  for  the  Ith  town  in  the  Jth  month. 

GWUSLD    -  The  amount  of  water  pumped  by  a  well. 

IDRZEN(l)  -  A  code  indicating  the  node  form  which  diversions  are  made  to 
a  reservoir  in  the  Ith  subbasin. 

IFILST(I)  -  A  code  indicating  the  subbasin  of  a  reservoir  to  which  water  is 
transferred  from  a  reservoir  in  the  Ith  subbasin. 

I FLAG     -  A  code  indicating  the  number  of  towns  that  could  not  meet  their 
demand  with  their  normal  sources. 

IFTIML(I)  -  The  number  of  months  of  recharge  that  will  be  computed  in  the 
Ith  month. 

IITIME    -  The  monthly  counter  ranging  from  I  to  the  total  number  of  months 
of  hydrologic  data. 

1MPMAX(I)  -  The  pump  capacity  at  an  imported  augmentation  source  in  the 
Ith  month. 

IMPMIH(I)  -  The  minimum  flow  that  must  be  exceeded  for  diversion  from  an 
imported  augmentation  source  in  the  Ith  month. 

IMPNOD    -  The  node  where  an  imported  augmentation  diversion  enters  the 
river. 

IMPREC    -  A  code  indicating  which  of  the  discharge  records  is  the  imported 
augmentation  source. 

Ih  -   The  standard  input  file  code. 

IORDER(I)  -  The  code  for  the  subbasin  of  the  reservoir  that  diverts  water 
from  the  river  in  the  Ith  order. 

IPRELM    -  The  number  of  months  in  an  initial  simulation  period  that 
stabilizes  reservois  and  well  effects. 

IRTURN(I)  -  The  node  to  which  the  sewage  of  the  Ith  town  returns  to  the 
streamflow. 

ITIME     -  The  monthly  counter  cycling  from  I  to  12. 
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ITITLE   -  Two  definitions: 

1)  In  the  mainline,  an  80  column  title  card  used  to  give  a 
brief  description  of  the  alternative  being  analyzed. 

2)  In  the  plot  routine,  a  descriptor  to  define  what  is  to  be 
plotted: 

to  plot  a  storage  level  hydrograph;  ITITLE  -  STOR 
to  plot  a  discharge  hydograph;  ITITLE  =  FLOW. 

ITRZFL   -  A  code  indicating  whether  a  discharge  trace  should  be  stored. 

ITRZST   -  A  code  indicating  whether  a  storage  level  trace  should  be  stored. 

IVAR(I)  -  A  code  indicating  which  of  the  four  time  series  is  to  be  plotted. 

LABEL (I)  -  The  name  of  the  Ith  variable  in  the  input  data  check. 

LINE (I)   -  The  Ith  position  on  one  line  in  the  hydrograph. 

MAXCAP(I)-  The  usable  storage  capacity  of  the  reservoir  at  the  1th  subbasin. 

MAXMDC(I)-  The  maximum  amount  of  water  obtained  by  the  Ith  town  from  some 
external  source  in  one  month. 

MAXST(I)  -  The  highest  storage  level  reached  by  the  reservoir  in  the 
Ith  subbasin. 

MDCCAP(I)-  The  capacity  of  the  Ith  town's  connection  to  some  external  source 

MDCUSE   -  The  amount  of  water  obtained  from  some  external  source. 

MINST(I)  -  The  lowest  storage  level  reached  by  the  reservoir  in  the 
Ith  subbasin. 

MONTH (I)  -  The  three  letter  abbreviation  for  the  Ith  month. 

NNODE    -  The  number  of  nodes. 

NOBS     -  The  number  of  observation  nodes. 

NODE     -  The  number  of  the  node  to  be  used  for  plotting  a  hydrograph. 

NODNAM(I)-  A  twelve  character  name  for  the  Ith  node. 

NPRINT   -  A  code  indicating  which  output  options  are  in  effect: 

if  evenly  divisible  by  2,  then  print  each  river  diversion  made; 

if  evenly  divisible  by  3,  then  print  each  reservoir  transer  made; 
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if  evenly  divisible  by  5,  then  print  each  imported  augmentation 
diversion; 

if  evenly  divisible  by  7,  then  print  a  list  of  the  input  data. 

NRCRDS   -  The  number  of  discharge  record  stations  used. 

OTIME    -  The  total  number  of  months  of  hydrologic  data. 

NTOWN    -  The  number  of  towns. 

NWELL    -  The  number  of  wells. 

OBS(I)   -  A  list  of  nodes  used  as  observation  nodes. 

OBSMAX(I)-  The  minimum  discharge  at  the  Ith  observation  node. 

OUT      -  The  standard  output  file  code. 

PERCNT(I,J)-  The  percentage  of  the  incremental  discharge  of  the  nearest 

downstream  gage  contributed  by  the  Ith  subbasin  in  the  Jth  month. 

PRIORY (I, J)-  A  code  listing  the  Jth  ordered  supply  source  for  the  Ith  town. 
The  list  is  In  two  groups,  the  first  eigt  being  normal  sources, 
the  second  eight  being  extended  sources. 

For  the  first  eight: 

1  indicates  wells; 

2  indicates  a  surface  water  diversion; 

any  number  greater  than  2  indicates  the  desired 
reservoir's  subbasin  number  plus  10. 

For  the  second  eight: 

1  indicates  an  external  source; 

2  indicates  wells; 

any  number  greater  than  2  indicates  the  desired 
reservoir's  suboasin  number  plus  10. 

Q(I)  -  The  discharge  at  the  Ith  node. 

QIMP  -  The  amount  of  the  imported  augmentation  discharge. 

QQ(I)  -  The  amount  of  the  imported  augmentation  discharge  at  the  Ith  node. 

RAInjGE  -  The  maximum  value  of  the  time  series  to  be  plotted. 

RATE  -  The  percentage  the  demand  is  of  a  town's  well  capacity. 
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RECORD (I)  -  The  discharge  record  that  the  discharge  of  the  Ith  node  is 
computed  from. 

RESDA(I)   -  The  drainage  area  of  the  reservoir  in  the  Ith  subbasin. 

RESOOB(I)   -  A  code  indicating  whether  the  reservoir  of  the  Ith  subbasin 
is  outside  the  basin's  drainage  boundary. 

RESQ       -  The  discharge  that  flows  into  a  reservoir . 

'RTURN(I,J)  -  The  percentage  of  water  used  by  the  Ith  town  that  returns 
to  streamflow  in  the  Jth  month. 

RUNOFF     -  The  average  of  the  monthly  runoff  (mgd/sq  mi)  of  each  discharge 
record  station  read  in. 

SCALE      -  The  scaling  interval  printed  on  the  plot. 

STCAPY(I)   -  The  storage  level  of  the  reservoir  in  the  Ith  subbasin. 

STMAX(I,J)  -  The  pump  capacity  for  removing  water  from  the  reservoir  in 
the  Ith  subbasin  in  the  Jth  month. 

STUSED     -  The  amount  of  water  taken  from  reservoir  storage. 

SUBADA(l)   -  The  drainage  area  of  the  Ith  subbasin. 

SUPPLY     -  The  amount  of  water  taken  from  a  supply  system  toward  meeting 
a  town's  demand. 

SWMAXCljJ)  -  The  pump  capacity  for  diverting  water  directly  to  the  Ith 
town  without  storage,  in  the  Jth  month. 

SWMIN(I,J)  -  The  minimum  flow  that  must  be  exceeded  before  diversion  at  the 
Ith  node  in  the  Jth  month. 

SWMINF     -  A  unit  discharge  (mgd/sq  mi)  used  to  compute  minimum  flows. 

SWSMAX(I,J)  -  The  pump  capacity  for  diverting  water  to  the  reservoir  in 
the  Ith  subbasin  in  the  Jth  month. 

SWUSED     -  The  amount  of  surface  water  diverted  for  use  without  storage. 

TOWN (I)    -  The  node  number  of  the  Ith  town. 

TTLDF(I)   -  The  total  deficits  accumulated  by  the  Ith  town. 

TTLDV(I)   -  The  total  diversions  accumulated  by  the  reservoir  in  the 
Ith  subbasin. 
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TTLMDC(I) 


-  The  total  amount  of  water  obtained  form  some  external  source 
by  the  Ith  town. 


TYMFAC      -  A  scale  factor  listing  a  conversion  for  the  rate 
(30.44  x  mgd  «=  mg  per  month) 

UNIT (I, J)   -  A  unit  function  that  account  for  the  affect  on  a  stream 
caused  by  pumping  the  Ith  well  in  the  Jth  future  month. 

WELCAP(I,J,K)  -  the  pump  capacity  of  the  Jth  well  owned  by  the  Ith  town 
in  the  Kth  month. 


WELL (I, J) 
WELNOD(I) 
WELPCT(I) 

X(l) 

XCHANG 

Y(I) 


-  the  well  number  of  the  Jth  well  owned  by  the  Ith  town. 

-  The  subbasin  number  in  which  the  Ith  well  is  located. 

-  The  percentage  of  the  subbasin 's  drainage  area  upstream  from 
the  Ith  well. 

-  The  Ith  preset  range  interval. 

-  The  amount  of  water  transferred  from  one  reservoir  to  another 

-  The  Ith  preset  scale  factor  associated  with  the  respective 
X  range  interval. 
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APPENDIX  II 
PROGRAM 

The  following  is  a  listing  of  the  Fortran  IV  source  deck  of  the  model. 


C  IPSWICH  SIMULATION  MODEL  -  VERSION  IB   FIXED  AUGUST  1971 

C 

C  SEE  U.S. GEOLOGICAL  SURVEY  OPEN  FILE  REPORT  'COMPUTER 

C  SIMULATION  MODEL  OF  THE  IPSWICH  RIVER  BASIN',  BURNS, A. W.  AND 

C  JAMES, I.C. II  FOR  EXPLANATION  OF  VARIABLES  AND  METHODOLOGY  OF 

C  THE  MODEL. 

C 

C  THIS  PROGRAM  WAS  CODED  BY  A.W.  BURNS,  U. S. GEOLOG ICAL  SURVEY, 

C         WATER  RESOURCES  DIVISION,  BOSTON,  MASSACHUSETTS. 

C 

C 

C 

C  THE  MAINLINE  MERELY  INITIALIZES  VARIABLES  AND  CALLS  THE 

C  SUBROUTINES. 

C 

C 

CCMMON     I   •  ITIME,  MAXCAP(60),  NNODE,  Q(60),  STCAPYI60), 
1  TOWNC25),  TYMFAC 

COMMON  /M12/  DA(IO),  IMPMAXU2),  IMPMINU2),  IMPNOD,  IMPQ,  IMPREC, 

1  IN,  NRCRDS,  NWELL,  PERCNTI60  ,  12 ) ,  00(60),  QSUBAI60),  REC0RDC60), 

2  RESDAI60),  RES008(60),  SUBACA(60),  WELN0DI60),  WELPCTI60) 
COMMON  /M35/  DEFICT<25),  DEMAND! 25 , 1 2  I ,  GWMAXI 25 , 12  ) ,  IFLAG, 

1  IPRELM,  IRTURN<25),  M0CCAPI25I,  NTIME,  NTOWN,  PR IORY( 25, 16) , 

2  RTURN(25,12),  SWMAX ( 25  ,12 1 ,  WELCAPC 25 , 10 ♦ 12 ) 

COMMON  /M78/  CLASS(7),  FRE0(10,8),  MAXSTI60),  MINST(60),  NOBS, 
1  OBS(IO),  OBSMAX(IO),  OBSMIN(IO) 

COMMON  /M568/  IDVRZN(601,  IFILST(60),  MAXMDCI25),  NODNAM (63,3 ) , 
1  OUT,  SWSMAX(60,12),  TTLDF(25I,  TTLDV(60),  TTLMDC(25),  I0RDER(60) 

COMMON  /M24/  EFF(60,12),  UNIT(60,12),  WELL(25,10) 

REAL  IMPMAX,  IMPMIN,  IMPQ 

INTEGER  RECORD,  RESOOB,  WELNOD 

REAL  MDCCAP 
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Appendix  2  -  Fortran  IV  Source  Deck 


C 
C 

C 
C 

c 

c 
c 


c 
c 

c 

c 
c 
c 


c 

c 
c 


c 
c 
c 


c 
c 
c 
c 


INTEGER  PRIORY 
REAL  KAXST,  MINST 
INTEGER  ORS 
REAL  MAXMDC 
INTEGER  CUT 
REAL  PAXCAP 
INTEGER  TOWN 
INTEGER  WELL 
DIMENSION  ITITLE(20» 


READ  (5,1000)  IN,  OUT 

READ  THE  TITLE  OF  THE  DATA  AND  THEN  PRINT  IT 

READ  (IN, 1010)  ITITLE 
INITIALIZATION 

DO  80  J  =  1,60 
80  ICPOER(J)  =  0 

CALL  SUBROUTINE  ■RDPARM*  TC  READ  ALL  THE  PARAMETERS 

CALL  RDPARM 

INITIALIZATION 

DO  60  L  =  l.NNODE 
OC(L)  =  0. 
DO  60  J  =  1,12 
60  EFF(L,J»  =  0. 

RUN  AN  INITIAL  PERIOC  TO  STABILIZE  RESERVOIRS  AND  WELL  EFFECTS 

DC  70  I  =  1, IPRELM 
ITIME  =  I  -  (( I-l)/12)*12 
CALL  RDFLOW 
CALL  NATHYO 
CALL  LOCAL 
70  CALL  FILLST 

INITIALIZATION 

DO  20  L  =  I, NOBS 

OBSMAX(LI  =  0. 

OBSMIML)  =  10000. 

DO  20  J  =  1,8 
20  FRE0(L,J1  =  0. 

DO  30  L  =  1,NT0WN 

TTLDF(L»  =  0. 

TTLMDC(L)  =  0. 
30  MAXMDC(L)  =  0. 

DO  40  L  =  l.NNODE 

TTLDV(L)  =  0 

MAXST(L)  =  0. 
40  MINST(L)  =  10000. 

INCREMENT  THE  MONTHLY  TIME  INTERVAL  THROUGH  THE  ENTIRE  RECORD 
OF  HYDROLOGIC  DATA 

DO  10  I  =  1,NTIME 

ITIME  =  I  -  (( I-l)/12)*12 

CALL  RDFLOW 
CALL  KATHYD 
IFLAG  =  0 
CALL  LOCAL 
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IF  (IFLAG.NE.OI 
CALL  FILLST 
CALL  KFPTRK 
CALL  APRINT 

STOP 


CALL  REGICN 


FORMAT 
FORMAT 
FORMAT 
END 


(211) 

(20AM 

(1H1  ////25X, 20A4////) 


C 

c 

C 

C 

c 
c 


c 
c 


SUBROUTINE  RCPARM 


READS  ALL  OF  THF 
IN  ORDER 


CCMMON  IITIMFt  ITIME, 
1  TOWN(25),  TYM<=AC 
CCMMON  /M12/  DAllClt  IMPMAXI12), 

1  IN,  NRCPDS,  NWFLL,  °FRCNT ( 60 , 1 2 ) 

2  RFSDM60)  ,  RES0GB(60)  , 
COMMON  /M?5/  OEFICK25) 


PARAMETERS  ANC  CHECKS  TO  SEE  THAT  THEY  APE 


MAXCAP(60),  NNODE,  Q(60),  STCAPY(6C) 


IVPVTMH2),  IMPNOD,  IMPG,  IMDREC, 
,  0Q(60)t  QSUBA(60),  PFC0Rn(6n), 
SURACA(60),  WELN0C(60),  WELPCT(60) 
DEMAND(25,12)  ,  GWMA X ( 25 , 1 2 ) ,  IFLAG, 


PRIOPYl 25,16) t 

NOBS, 


1  IPRPLM,  IRTURN125),  MDCCAPJ25),  NTIME,  NTOWN, 

2  RTURN( 25,12) ,  SW^AX ( 25  , 1 2 ) ,  WELC AP ( 25 , 10 , 12 ) 
CCMMON  /M78/  CLASSI7),  FRE0(10t8)«  MAXST(60),  MINST(60) 

1  CBS(IO),  CBSMAX(IO),  CESMIN(IO) 

COMMON  /M568/  I0VRZN(60),  IFILSTI60),  MAXMDC(25),  NODNAM ( 63 , 3 ) , 
1  OUT,  SWSMAX(60, 12),  TTLDF(25),  TTLDV(60),  TTLM0C(25),  ICRDER(60) 
CCMMON  /S2356/  NPPINT,  ST  MAX ( 60, 12 ) ,  SWMIN(60,12) 
COMMON  /M24/  EFF(60,12),  UNIT(60,12),  WELL(25,10) 

ITRZFL,  ITRZST,  CTRACE ( 1 20 , 9 ) ,  FLTR AZ (  1  20  ,9 )  ,  IWANTX 
IMPMIN,  IMPQ 


CCMMCN  /R2/ 
REAL  IMPMAX 

MDCCAP 

PAXST,  MINST 

MAXMOC 

MAXCAP 

RECORD, 


REAL 

Ri=AL 

RFAL 

REAL 

INTEGER 

INTEGER 

INTEGER 

INTEGER 

INTEGER 

INTEGER 


WELNOC 


RESOCP, 
PRIORY 
CBS 
CUT 
WELL 

TCWN,  TITLE 
DIMENSION  LABELI32),  CCNSUP(50) 

DATA  LABEL  /  4H     , 4HI MPC, 4HREC0,4HWELN ,4HT0WN, 4H0R S  ,4HIRTU, 
14HPERr,4HRESC,4HSURA,4HIFIL,4HMAXC,4FSWMI,4HSWMA,4HSWSM,4HSTMA, 
24HRESC,4HDEMA,4HRTUR,4HMCCC,4HPRI0,4HWELC ,4HWELP ,4HWELL , 4HI MPM, 
34HIMPM,4HCLAS,4HTYMF,4HDA   , 4HCCNS ,4HUNI T ,4H I DVR  / 


L  =  1 
LL  =  1 

READUN,1010)TITLE,  NUV,  KCDE,  NNCDE,  NTIME, 
1  NRCRDS, IPRELM,NPRINT, ITRZFL, ITRZST,  IWANTX 
READ  (IN, 1040)  ( (N0DNAM(I,J),J=1,3), 1=1, NNODE) 
IF  (TITLE. NE. L'ABEL(L)  )  GO  TO  700 
L  =  2 

READ  (IN, 1010)  TITLF,NUM, IMPREC, IMPNOD 
IF  (TITLE. NE. LABEL(L) )  GC  TO  700 
L  =  3 
LL  =  0 

DC  10  M  =  1, NNCDE, 24 
LL  =  LL  +  1 
L  =  3 

LL  =  0 
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DC  10  M  =  1,NN0DE,24 

LL  =  LL  ♦  1 

N  =  M  ♦  23 

IF  (N.GT.NNODE)  N  =  NNOCE 

READ  (IN, 1010)  TITLE, NUM  ,  ( RECORD! 1 1 » I  =  M,N) 

IF  (TITLE. NE.LABEL(L) .OR. NUM. NE.LL)  GO  TO  700 
10  CONTINUE 

L  =  4 

LL  =  0 

DC  20  M  =  1,NWELL,24 

LL  =  LL  +  1 

N  =  M  +  23 

IF  (N.GT.NWELL)  N  =  NWELL 

READ  (IN, 1010)  TITLE, NUM,  ( WELNOD( I ) , 1=  M,N) 

IF  (TITLE. NE.LABEL(L). OR. NUM. NE.LL)  GO  TO  700 
20  CONTINUE 

L  =  5 

LL  =  C 

DO  30  f>    =  1,NT0WN,24 

LL  =  LL  +  1 

N  =  M  +  23 

IF  (N.GT.NTOWN)  N  =  NTOWN 

READ  (IN, 1010)  TITLE, NUM,  (TOWN ( I ) , I  =  M, N ) 

IF  (TITLE. NE.LABEL(L).CP. NUM. NE.LL)  GO  TO  700 
30  CONTINUE 

L  -    6 

LL  =  1 

READ     (IN, 1010)    TITLE, NUM.     ( 0BS( I ) , I = 1 ,NOBS) 

IF    (TITLE. NE. LABEL(LI )    GO    TO    700 

L    =    7 

LL    =    0 

DO  40  P  -  1,NTCWN,24 

LL  =  LL  ♦  1 

N  =  M  ♦  23 

IF  (N.GT.NTOWN)  N  =  NTOWN 

READ  (IN, 1010)  TITLE, NUM,  (IRTURN(I),  I=M,N) 

IF  (TITLE. NE.LABEL(L). OR. NUM. NE.LL)  GO  TO  700 
40  CONTINUE 

L  =  8 

DC  50  LL  =  1,NN0DE 

READ  (IN, 1020)  TITLE,  NU*,  ( PERCNT ( LL, J) , J  =  1,12) 

IF  (TITLE. NE.LABEL(L). OR. NUM. NE.LL)  GO  TO  700 
50  CCNTINUE 

L  =  9 

LL  =  0 

OC  60  M  =1,NN0DE,12 

LL  =  LL  +  1 

N  =  M  ♦  11 

IF  (N.GT.NNODE)  N  =  NNOCE 

READ  (IN, 1030)  TITLE, NUP.  (RESDA(I),I  =  M,N) 

IF  (TITLE. NE.LABEL(L). OR. NUM. NE.LL)  GO  TO  700 
60  CONTINUE 

L  =  10 

LL  =  0 

OC  70  M  =  1,NN0DE,12 

LL  =  LL  +  1 

N  =  M  ♦  11 

IF  (N.GT.NNODE)  N  =  NNOCE 

READ  (IN, 1030)  TITLE,  NUM,  (SUBADA(I),I  =  M,N) 

IF  (TITLE. NE.LABEL(L). OP. NUM. NE.LL)  GO  TO  700 
70  CONTINUE 

L  =  11 

LL  =  0 

DC  80  M  =  1,NNCDE,24 

LL  =  LL  +  1 

N  =  M  +  23 

IF  (N.GT.NNODE)  N  *  NNOCE 

READ  (IN, 1010)  TITLE, NUM, ( IF  I LST( I ) , I=M,N ) 
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IF  (TITLE. NE.LABEL(L). OR. NUM. NE.LL)  GO  TO  700 
80  CONTINUE 

L  =  12 

LL  =  0 

OC  <30  M  =  1,NNCDE,12 

LL  =  LL  ♦  1 

N  =  f    +  11 

IF  (N.GT.NNODE)  N  =  NNOCE 

READ  (IN, 1030)  TITLE,  NL'M,  (MAXCAP(I),I  =  M,N) 

IF  (TITLE. NE. LABEL(L) .CP. NUM. NE.LL)  GO  TO  700 
90  CONTINUE 

DC  95  I  =  l.NNODE 
95  STCAPY(I)  =  ,5*MAXCAP(I) 

L  =  13 

LL  =  1 

READ  (IN, 1020)  TITLE, NUM,  SWMINF 

IF  (TITLE. NE.LABEL(LI)  GO  TO  700 

AREA  =  0 

DC  IOC  I  =  1,NN0DE 

AREA  =  AREA  ♦  SUBA0A(  I) 

X  =  SWMINF*AREA 

OC  100  J  =  1,12 
100  SWMIN(I,J)  =  X 

L  =  14 

DO  11C  LL  =  l.NTOWN 

READ  (IN, 1030)  TITLE,  NUM.  ( SWMAX ( LL, J) , J  =  1,121 

IF  (TITLE. NE.LABELCL). OP. NUM. NE.LLI  GO  TO  700 
110  CONTINUE 

L  =  15 

DO  120  LL  =  1,NNC0E 

READ  (IN, 1030)  TITLE,  NUM,  ( SWSMAX (LL . J ) , J  =  1,12) 

IF  (TITLE. NE.LABEL(L). OR. NUM. NE.LLI  GO  TO  700 
120  CCNTINUE 

L  =  U 

DC  130  LL  =  l.NNODE 

READ  (IN, 1030)  TITLE,  NUM,  (STMAX( LL , J  I , J  =  1,12) 

IF  (TITLE. NE.LABEL(L). OR. NUM. NE.LLI  GO  TO  700 
130  CCNTINUE 

L  =  17 

LL  =  0 

DC  140  M  =    l.NNODE, 24 

LL  =  LL  ♦  1 

N  =  M  +  23 

IF  (N.GT.NNODE)  N  =  NNODE 

READ  (IN.IOIOI  TITLE, NUM,  ( RE SOOBf 1 1 , I =M ,N) 

IF  (TITLE. NE.LABEL(L). OR. NUM. NE.LL)  GO  TO  700 
140  CONTINUE 

L  =  18 

DC  150  LL  =  l.NTCWN 

READ  CIN, 10201  TITLE, NUM,  ( DEMAND! LL, J ) ,  J  =  1,121 

IF  (TITLE. NE.LABEL(L). OR. NUM. NE.LLI  GO  TO  700 
150  CCNTINUE 

L  =  IS 

DC  160  LL  =  l.NTOWN 

READ  (IN, 1020)  TITLE, NUM,  (  RTURN(LL,J),  J=  1,12) 

IF  (TITLE. NE.LABEL(L) .OR. NUM. NE.LL)  GO  TO  700 
160  CCNTINUE 

L  =  20 

LL  =  0 

00  170  M  =  l.NTOWN, 18 

LL  =  LL  ♦  1 

N  =  M  ♦  17 

IF  (N.GT.NTOWNI  N  =  NTOWN 

READ  (IN, 1020)  TITLE, NUK,  (MOCCAP(I),I  =  M,N) 

IF  (TITLE. NE.LABEL(L). OR. NUM. NE.LL)  GO  TO  700 
170  CONTINUE 

L  =  21 

DO  180  LL  =  1,  NTOWN 
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READ  (IN, 1010)  TITLE, NUN,  ( PR IORY ( LL , J  )  . J  =  1,16) 

IF  (TITLE. NE.LABEL(L) .OR. NUM. NE.LL)  GC  TO  700 
180  CCNTINUE 

L  =  22 

LL  =  0 

DC  195  I  =  1,NT0WN 

DC  190  J  =  1,10 

LL  =  LL  +  1 

IF  (LL.GE.100)  LL  =  0 

READ  UN.1C20)  TITLE,  Ml*,  I WELCAPU  ,  J,K)  ,  K  =  1,12) 

IF  (TITLE. NE.LABEL(L). OR. NUM. NE.LL)  GO  TO  700 

IF  (WELCAP(I,J,l).EQ.O. )  GO  TO  192 
190  CCNTINUE 

GO  TO  195 

192  OC  193  KK  =  J, 10 
DC  193  K  =  1,12 

193  WELCAPI I,KK,K)  =  0. 
195  CCNTINUE 

L  =  23 

LL  =  C 

DC  200  M  =  1, NWELL, 18 

LL  =  LL  +  1 

N  =  f*    +  17 

IF  (N.GT. NWELL)  N  ■  NWELL 

READ  (IN,1C20)  TITLE, NUM.  (WELPCT(I1,I  =  f,N» 

IF  (TITLE. NE.LABEL(L». OP. NUM. NE.LL)  GO  TO  700 
200  CONTINUE 

L  =  24 

DC  210  LL  =  l.NTOWN 

READ  (IN, 1010)  TITLE, NUM,  < WE LL ( LL, J  I , J  =  1,101 

IF  (TITLE. NE.LABEL(L) .OR. NUM. NE.LL)  GO  TO  700 
210  CCNTINUE 

L  =  25 

LL  =  1 

READ  (IN, 1030)  TITLE, NUP,  UMPMAXU),  I  =  1,12) 

IF  (TITLE. NE. LABEL(L) )  GO  TO  700 

L  =  26 

READ  (IN, 1030)  TITLE, NUP,  (IHPHIN(I),  I  =  1,12) 

IF  (TITLE. NE. LABELCLI )  GO  TO  700 

L  =  27 

READ  (IN,1C30)  TITLE, NUK,  (CLASSCII,  I  =  1,7) 

IF  (TITLE. NE. LABEL(L) )  GO  TO  700 

L  =  28 

READ  (IN, 1030)  TITLE, NUK,  TYMFAC 

IF  (TITLE. NE. LABEL(L) )  GO  TO  700 

L  =  29 

READ  (IN, 1030)  TITLE, NUH,  (DA(I),I  =  l,NRCRDS) 

IF  (TITLE. NE.LABEL(D)  GO  TO  700 

DO  220  I  =  1,NTCWN 

DC  22C  K  =  1,12 

GWMAXU.K)  =  0. 

DC  22C  J  =  1,10 
220  GWMAXU.K!  =  GWMAXU.K)  ♦  WELCAPU,J,K) 

L  =  30 

LL  =  C 

00  230  M  =  l,NT0WN,18 

LL  =  LL  ♦  1 

N  =  ¥    +  17 

IF  (N.GT.NTOWN)  N  =  NTCWN 

READ  (IN, 1020)  TITLE,  NLM,  (CONSUPf I ) , I=M,N) 

IF  (TITLE. NE.LABEHLI.CR. NUM. NE.LL)  GO  TO  700 
230  CCNTINUE 

L  =  31 

DC  235  LL  =  1, NWELL 

READ  UN, 1030)  TITLE,  NUM,  ( UNIT! LL , I ) , I  =  1,12) 

IF  UITLE.NE.LABEL(L).CR.NUM.NE.LL)  GO  TO  700 
235  CONTINUE 

L  =  32 
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DO  236  M  =  1,NN00E,24 
LL  =  LL  ♦  1 
N  =  M  *    23 

IF  (N.GT.NNODE)  N  =  NNOCE 

READ     (IN, 1010)    TITLE,     NUM,     I IOVRZNf 1 1 1 1    =    M,N) 
IF     (TITLE. NE.LABEL(L). OR.    NUM.NE.LL)    GO    TO    700 
?36    CCNTIhUE 
J    =    0 

237  I H I K    =    NNGCE 

DC    238    I    =    l.NNODE 

IF    ( IDVRZNU J.GT.O.AND. ICVRZN(I).LT. IMIN)     IMIN    =     I 

238  CONTINUE 

IF  ( IPIN.EC.NNODEI  GC  TC  239 
J  =  J  ♦  1 
ICRCER(J)  =  IMIN 
IDVRZMIMIN)  =  -IOVRZM  IMIN) 
GO  TO  237 

239  DC  241  I  =  1,NNCDE 

241  IOVRZMI)  =-IDVRZN(I) 
OC  240  I  =  l.NTOWN 
DO  240  J  =  1,12 

240  DEMANO(I,J»  =  CONSUP (  I ) *DEMANC ( I , J  I / . 83333 
C 

C         WITH  PROPER  CPTICN  CCDE  CALL  THE  ROUTINE  'LST»  TO  LIST  THE 

C  INPUT  OATA. 

C 

IF  ( (NPRINT/7I*7.EQ.NPRINT)  CALL  LST 
RETURN 
700  WRITE  (OUT, 2000)  LL,LABEL(L) 
STCP 
C 
C 
1010  FCRVAT  (A4.2X, 12,2413) 
1020  FORMAT  ( A4 , 2X, I  2, i 8F4. 3  ) 
1030  FORMAT  ( A4 , 2X , I  2, 12F6  .  1  ) 
1040  FORMAT  (3A4) 

2000  FORMAT  (10X,«***  CARD  ,,I3,«  CF  ■ , A4 , •  IS  OUT  OF  CRDER  •) 
END 


SUBROUTINE  RCFLOW 
C 
C 

C  DISTRIBUTES  THE  STREAMFLOW  TO  EACH  SUBBASIN 

C 
C 

COMMON  IITIME,  ITIME,  MAXCAP(60),  NNCOE,  C(60),  STCAPY(60), 
1  T0WN(25),  TYMFAC 

COMMON  /M12/  DA(IO),  IMFMAXU2),  IMPMINI12),  IMPNOD,  IMPQ,  IMPREC, 

1  IN,  NRCRDS,  NWELL,  PERCNT(60 ,12) •  00160).  QSUBA160),  REC0RDI60), 

2  RESDA(60),  RES0CB(60),  SUBADA(60>,  WELN0DC60),  WELPCTI60) 
REAL  IMPMAX,  IMPMIN,  IMPQ 

INTEGER  RECORD,  RESOOB,  WELNOO 
REAL  MAXCAP 
INTEGER  TCWN 
DIMENSION  DIS( 10) 
C 

READ  (IN, 1000)  (OISC I), 1  =  1,  NRCRDS) 
C 
C         COMPUTE  THE  AVERAGE  RUNOFF  FOR  FOP  ALL  THE  GAGING  STATIONS  USED 


C 


RUNOFF  =  0 
DO  10  I  =  1, NRCRDS 
10  RUNOFF  =  RUNOFF  ♦  DIS(I)/0A(I) 

RUNOFF  =  RUNOFF*. 646317/  FLOAT( NRCRDS  ) 
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C  DISTRIBUTE  THE  INCREMENTAL  GAGED  STREAMFLCW  TO  EACH  SUBBASIN 

C  ANC  PUT  AVERAGE  RUNOFF  INTO  EACH  RESERVOIR. 

C 

K  =  1 

DC  20  INOCE  =  l.NNODE 
I  =  RECCRD(INCOE) 
IF( I .EQ.K)  GO  TO  15 
K  =  I 
L  =  I  -  1 
12  DIS(  II  =  DISU)  -  DIS(L> 
L  =  L  -  1 

IF  (L.NE.O)  GC  TO  12 
15  OSUBA(INCDE)  =  PERCNT ( I NOCE, I TI WE  I *D IS ( I  )  *. 64631 7 
PESQ  =  0 

IF  (SUBADA(INODE).NE.O.  1  PESC  =  (RESCA{ INODE l/SUBADA( INODE ))*QSUBA 
1UN0DE) 
IF  (RESOOB(INODE).NE.l)  QSUBA(INODE)  =  QSUBA(INODE)  -  RESQ 
IF  (RESOCB(INCDE).EQ.ll  RESG  =  PUNCFF*RESDA( INODE ) 

STCAPYI INODE)  =  STC APY ( I  NODE)  ♦  RESC*TYMFAC 
IF  (STCAPY(IN0DE).LE.MAXC/IP(  INODE)!  GO  TO  20 

IF  (RESOOB(INODE).NE.l)  QSUBA(INODE)  =  QSUBAIINODE)  4  (STCAPYC INOD 
IE)  -  MAXCAPi INODE) )/  TYMFAC 
STCAPYt INCDE)  =  MAXCAP (  INCDE ) 
20  CCNTINUE 

IMPG  =  DISUMPREO*. 646317 
RETURN 
C 
C 
1000  FORMAT  (  1CF8.0) 
END 


SUBROUTINE  NATHYC 
C 
C 

C         COMPUTES  THE  EFFECTS  OF  WELLS  ON  THE  STREAMFLOW  AND 
C  ACCUMULATES  THE  STREAMFLOW  OF  EACH  SUBBASIN  IN  THE  RIVER. 

C 
C 

CCMMCN  IITIME,  ITIME,  MAXCAP(60>,  NNODE,  Q(60),  STCAPY160), 
1  T0WN(25),  TYMFAC 

COMMON  /M12/  DA(IO),  IMPMAXU2),  IMPMINI12),  IMPNCO,  IMPQ,  IMPREC, 

1  IN,  NRCRDS,  NWELL,  PERCNT ( 60 , 1 2 ) »  00(601,  QSUBAI60),  REC0RD(60), 

2  RESDAI60),  RES00BC60),  SUBADA(60),  WELN0DI60),  WELPCT(60) 
COMMON  /S2356/  NPRINT,  STMAXI 60, 12 ) ,  SWMIN(60,12) 

CCMMON  /M24/  EFF(60,12),  UNITI60.12),  WELL(25,10) 

COMMON  /M568/  IDVRZN(60),  IFILSTI60),  MAXMDCC25),  NODNAM (63, 3 ) , 
1  OUT 

INTEGER  CUT 

REAL  IMPMAX,  IMPMIN,  IMPQ 

INTEGER  RECORD,  RESOOe,  WELNOD 

INTEGER  WELL 

REAL  MAXCAP 

INTEGER  TOWN 

DIMENSION  IFTIME14) 

DATA  IFTIME  /0,2,6,12  / 
C 
C 

C         COMPUTE  THE  AMOUNT  OF  AN  IMPORTED  DISCHARGE  AND  AUGMENT  THE 
C  STREAMFLOW  THAT  AMOUNT   AT  THE  PROPER  RIVER  POINT 

C 

01  MP  =AM  INK  IMPMAX  (  IT  IME) , IMPQ-IMPMIN(  IT  IME)  ) 

IF  (QIMP.LT.O.)  OIMP  *  0. 

QQ(IMPNOD)  =  GIMP 

OIMP  =  OIMP  *TYMFAC 
C 
C         WITH  PROPER  OPTION  CCOE  PRINT  THE  AMOUNT  OF  IMPORTED  WATER  THAT 

29 


1 
I 

1 
1 

I 
I 
I 
I 
I 
I 

I 

1 
I 
I 

I 

I 
fj 


C  AUGMENTS  THE  RIVER  FLOW  AppendiX  *   -  FortMn  "  Wfg  °eCk 

C 

IF  ( (NPRINT/5)*5.EQ.NPRINT.AND.QIMP.NE.0.I  WRITE  (OUT, 2000)  QIMP  , 
1  I ITI ME 
C 

C         FOP  WELLS  AFFECTING  A  TRIBUTARY,  SUBTRACT  FROM  A  SUBBASIN'S 
C  RUNOFF  THE  EFFECT  CF  WELL  PUMPAGES  ANC  THEN  ACCUMULATE 

C  THE  RIVER  DISCHARGE  UP  TO  RIVER  PCINT  M« 

C 

I  =  1 
IWELL  =  1 
10  DC  50  INCDE  =  I, NNODE 

IF  (WELNOD(IWELL).NE. INCDE)  GC  TO  40 
IF  (WELPCTUWELL).GT.l.  )  GO  TO  60 
QDELTA  =  0 

DELTAC  =  OSUBAI I NOOE) *WELPCT (I  WELL )  -  EFFUWELL, 1) 
20  IF  IDELTAO.LT. 0.1  EFFUWELL, 2)  =  EFFUWELL, 2)  -  DELTAQ 
DELTAQ  =AMAX1(DELTAQ,0.) 
GCELTA  =  CCELTA  ♦  DELTAQ 
IUELL  =  IWELL  ♦  1 

IF  (WELNCDUWELL).NE. INCDE)  GC  TO  30 

DELTAC  =  QSUBAUNODE)*(WELPCTUWELL)  -  WELPCU I WELL-1 >) 
1  -  EFF( IWELL, 1 ) 
GO  TO  20 
30  QSUBA(INODE)  =  OSUBAI INODE )*( 1.  -  WELPCTC IWELL- 1)  I  ♦  QDELTA 
40  IF  (INODE. GT.l)  GO  TO  50 
0(1)  =  QSUBAU) 
GC  TO  70 
50  Q(INODE)  =  0<INCDE-il  ♦  QSUBA(INODE)  ♦  QQUNODEI 
GO  TO  75 
C 

C  FOR  WELLS  AFFECTING  THE  MAINSTREAM,  ACCUMULATE  THE  RIVER  DISCHARGE 

C  UP  TO  RIVER  POINT  »I"  AND  THEN  SUBTRACT  FROM  THAT  STREAMFLOW 

C  THE  EFFECT  OF  WELL  PUMPAGES. 

C 

60  Q(INOCE)  =  QUNODE-1)  +  QSUBAUNODE)  ♦  QQ(INODE) 
DELTAC  =  QUNCDE)  -  EFFUWELL, 1) 

IF  (DELTAO.LT.O.)  EFFUWELL, 2)  =  EFFUWELL, 2)  -  DELTAQ 
Q(INOCE)  =AMAX1(DELTAQ,0.I 
IfcELL  =  IWELL  ♦  1 
70  I  =  INODE  ♦  1 

IF  (I  .LE.'NNOOE)  GO  TO  10 
C 

C         IN  THE  SPRING  MONTHS  REDUCE  THE  WELL  EFFECT  HYDROGRAPHS  TO  ZERC 
C  ANC  SUBTRACT  THAT  AMOUNT  FROM  THE  STREA«FLOW  AT  EACH 

C  RESPECTIVE  RIVER  POINT. 

C 

75  IF  UTIME.EQ.1.CR.ITIME.GT.4)  RETURN 
L  =  IFTIMEUTIMEI 
00  130  IWELL  =  l.NWELL 
INCDE  =  WELNCDI IWELL! 
IF  (WELPCTUWELLI.GT.l.)  GC  TC  100 
DC  90  I  =  1VL 

IF  UCSUPAUNODE)*WELPCT(  IWELD/2.  )  .LE.  EFFUWELL,  I ) )  GO  TO  130 
QSUBAUNODE)  =  QSUBAUNODE)  -  EFFUWELL, II 
DO  80  K  =  INODE, NNODE 
80  Q(KT  =  0(K)  -  EFFUWELL, I) 
90  EFFUWELL, I)  =  0 
GC  TO  130 
100  DO  12C  I  =  1,L 

IF  <(CUN00E>/2.  I.  LE.  EFFUWELL,  It)  GC  TO  130 
DO  110  K  =  INODE, NNODE 
110  Q(K>  =  Q(KI  -  EFFUWELL, I) 
120  EFFUWELL,!)  =  0 
130  CONTINUE 
RETURN 


C 
C 


2000  FORMAT  (30X,F7.1,»  MG  WAS  IMPORTED  INTO  THE  BASIN  IN  MONTH  NUMBER 
1M3I 

END  30 
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SUBROUTINE    LOCAL 

C 

c 

C         DRAWS  WATER  FROM  FIRST  SET  OF  SUPPLY  SOURCES. 

C 

c 

COMMON  IITIME,  ITIME,  MAXCAP(60),  NNCOEt  C(60),  STCAPY(60), 
1  TOWN(25>,  TYMFAC 
CCMKON  /N35/  DEFICT(25),  0EMAND( 25, 12 ) ,  GWMAX ( 25, 12  I ,  IFLAG, 

1  IPRELM,  IRTURN(25),  MDCCAP(25),  NTIME,  NTOWN,  PRIORY(25  ,16 > , 

2  RTURN(25,12),  SWMAX ( 25, 12 ) •  WELCAP( 25, 10, 12 ) 
CCMMON  /S2356/  NPPINT,  STKAX (60. 12 ) ,  SWMIN(60,12) 
REAL  MAXCAP 

INTFGER  TCWN 

REAL  MDCCAP 

INTEGER  PRIORY 
C 
C 

DO  ICO  ITOWN  =  1, NTOWN 

INCDE  *  TOWN< ITOWNI 

SUPPLY  =  0 
C 

C  GO  TO  THE  SOURCE  OF  THE  HIGHEST  PRIORITY. 

C 

DO  70  K  =  1,8 

L  =  PRIORY(ITOWN.K)  ♦  1 

IF  (L.LE.3)  GO  TC  (10,20.30),  L 
C 

C  TAKE  WATER  FROM  A  RESERVOIR 

r 

ISTORE  =  L  -  11 

STUSED  =  AMIN1(STMAX( ISTORE, ITIME I , STCAPY ( I  STORE  J /TYMFAC ,  DEMANDd 
1TCWN, ITIPE)  -  SUPPLYI 
STCAPY(ISTORE)  =  STC APY ( I  STORE )  -  STUSED*TYMF AC 
SUPPLY  =  SUPPLY  ♦  STUSEC 
GC  TC  60 
C 

C  NO  OTHER  SOURCES  FROM  THE  FIRST  SET,  SO  DEFINE  A 

C  TEMPORARY  DEFICIT. 

C 

10  DEFICT( ITCWN!  =  DEMANCl  ITOWN,  ITIME )  -  SUPPLY 
IFLAG  =  IFLAG  ♦  1 
GO  TO  80 
C 

C  TAKE  WATER  FRCH  WELLS 

C 

20  RATE  *  (CEMANDCITOWN, ITIME)  -  SUPPLY ) /GWMAX( I  TOWN, ITIME I 
IF  (RATE. GT. 1.1  RATE  =  1. 
DO  25  I  =  1,10 

GWUSEC  =  RATE*WELCAP(  ITCWN, I, ITIME) 
IF  (GWUSED.NE.O.)  CALL  EFFECT (GWUSED, I , ITOWN ) 
25  SUPPLY  =  SUPPLY  ♦  GWUSEC 
GO  TO  60 
C 

C  TAKE  WATER  FROM  A  STREJMFLOW  DIVERSION 

C 

30  SWUSED  =  AMIN1(SWMAX( ITOWN, ITIME) ,Q( INODEI  -  SWMI N( I  NODE , ITI ME) , 
1  OEMANDUTOWN,  ITIME)  -  SUPPLY) 
IF  (SWUSED. LT.O.)  SWUSED  =  0. 
DC  40  1=  INOCE,NNODE 
40  0(1)  =  0( II  -  SWUSED 

SUPPLY  =  SUPPLY  ♦  SWUSED 
60  IF  (OEMANDUTOWN, ITIME). GE. (SUPPLY  -.000 1*SUPPL Y) . AND. DEMAND! ITOWN 

i,ITIME).LE. (SUPPLY  ♦  ,0001+SUPPLY ) )  GO  TO  75 
70  CONTINUE 
75  DEFICT(ITOWN)  =  0. 
C 
C         RETURN  USED  WATER  TO  STREAMFLOW  ACCORDING  TO  A  SEWAGE  FUNCTION. 
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80 


90 
100 


ADD  =  RTLRN( ITCWN,ITIME J+SUPPLY 
J  =  IRTURNUTOWN) 
DC  90  I  =  J.NNCOE 
0(1)  =  0(1)  ♦  AOD 

CONTINUE 

RETURN 

END 


C 

c 
c 
c 
c 
c 


c 
c 


SUBROUTINE  EFFECT( GWUSEC,  I  ,  I  TCWN ) 


ACCUMULATES  THE  EFFECTS  OF  WELL  PUMPAGES  ACCORDING  TO  A  UNIT 
FUNCTION  FOR  EACH  WELL. 


CCMMCN  IITIME,  ITIME,  MAXCAPI60),  NNODE,  Q(60),  STCAPY(60I, 
1  TOWNI25),  TYPFAC 
CCMMCN  /M24/  EFF(60,12>,  UNIK60.12),  WELL(25,10» 
INTEGER  TCWN 
INTEGER  WELL 
REAL  MAXCAP 


IWELL  =  WELLUTOWN.il 
DC  10  J  =  If  11 
10  EFF(IWELLtJ)  =  EFF ( I WELL . J+l  )  ♦  UNIT < IWELLt J  1 +GWUSED 
RETURN 
END 


C 

c 

c 
c 
c 


c 
c 


c 
c 
c 


SU8C0UTINE  REGION 


DRAWS  WATER  FROM  A  SECCND  SET  OF  SUPPLY  SOURCES, 


COMMON'  IITIME,  ITIME,  MAXCAP(60>,  NNCDE,  C(60),  STCAPY(60), 
1  T0WN(25),  TYMFAC 
COMMON  /N35/  DEFICT(25),  CEMAND ( 25 , 1 2 ) ,  GWMAX ( 25, 12 ) ,  TFLAO, 

1  IPRELM,  IRTURN125),  MDCCAP(25),  NTIME,  NTOWN,  PP I ORY ( 2 5 , 16 ) , 

2  RTURN(25,12),  SWMAX (  25 ,  1  2 ) ,  WELC AP ( 25 , 10, 12 ) 

COMMON  /M569/  IDVRZN160),  IFILST(60),  MAXMDC(25),  NODNAM { 63, 3 ) , 
1  OUT,  SWSMAX(60, 12) ,  TTLCF(25),  TTLOV160),  TTLMDC(25),  irP0ER(60) 
CC^MCN  /S2356/  NPRINT,  STMAX ( 60 , 1 2 ) ,  SWMIN(60,12) 
REAL  MAXCAP,MCCUSE 
INTEGER  TCWN 
RE*L  NDCCAP 
INTEGER  PRIORY 
REAL  VAXMCC 
INTEGER  CUT 


DC  70  ITCWN  =  1, NTOWN 
IF  (DEFICT(ITCWN).EQ.O.) 
INODE  =  TOWN! ITOWN)  -  1 


GO  TC  70 


GO  TO  THE  SOURCE  OF  THE  HIGHEST  PRIORITY, 


32 


I 
1 
I 
I 
I 
I 
I 

I 

I 

i 
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|_    =    PPIPRY(  I  TOWN  ,K)     *    1 

IF     (L.LE.3)     GO    TO     (60,10,20),    L 

r 

C  TAKE  WATER  FR€M  A  RESERVOIR 

r 

TST^RE  =  L  -  11 

STUSEC    =     AMIN1(STVAXUSTCPE,ITI*E),STCAPY{ I S TORE ) /TY ^FA C, CF F  ICT (  IT 
ICWN)) 

STCAPY( ISTCPE)     =    STCAPYUSTORE)     -    STUSED*TYMFAC 

DEFICK ITOWN)     =    OFF  ICT (  I  TTWN )     -    STUSED 

SUPPLY    =    STUSEC 

GC    TO    4  0 
C 

r  .  TAKE  WATER  F^CM  SO^E  EXTFRNAL  SOURCE 

C 

10  ACCUSE  =AMIN1(MDCCAP( ITWM ,DEFICT( ITCWN) ) 

DEPICT! ITCWN)  =  DEFICT{ ITCWN)  -  VDCUSE 

SUPPLY  -  ACCUSE 

^DCUSE  =  ^DCUSE  *  TYMFAC 

TTLvnC( ITCWN)  -  TTL^HCI  ITCWN)  +  MCCUSE 

MAXr'OCt  ITCWN)  =AM&xi  (MAX^CC(  I  TCWN  )  ,  MDCUS  E  ) 

GO  TT  40 

r 

C  take  WATER  FPCV  WELLS 

MAXfCC(  ITCWN)  =AMAX1 ( MAXHCCI I TCWN ) .MOCUSE I 
GO  TC  40 
C 

C  TAKE  WATER  FROM  WELLS 

C 

20  RATE  =  DEPICT! ITCWN) /GWMAX < I  TOWN, ITI ME ) 
SUPPLY  =  0 

IF  (RATE.GT.l.)  RATE  =  1. 
DO  25  I  =  1,10 

GWUSED  =  RATE*WELCAP( ITCWN, I , ITIME I 
IF  (GWUSED. NE.O.)  CALL  EFFECT (GWUSED, I, ITOWN ) 
25  SUPPLY  =  SUPPLY  ♦  GWUSED 

DEFICT(ITCWN)  =  DEF  ICT(  ITOWN )  -  SUPPLY 
C 

C         RETURN  USED  WATER  TO  STREAMFLOW  ACCORDING  TO  A  SEWAGE  FUNCTION. 
C 

40  ADD  =  RTURMITCWN, ITIME  )*SUPPLY 
J  -  IRTURN(ITOWN) 
DC  45  I  =  J.NNCDE 
45  0(1!  =  Q( I)  ♦  AOD 

IF  (DEFICTI  I  TO  V»N  ).  GT. -.000  LAND.  DEPICT  (I  TOWN  ).LT.. 0001)  DEPICT!  I  TO 
1WM  =  0. 
IF  (DEPICT(ITOWN).EO.O. )  IFLAG  =  IPLAG  -  1 
IF  (  IFLAG. EQ.O)  RETURN 
IF(DEFICT< ITCWN). EO.O.)  GC  TO  70 
50  CONTINUE 
C 

C         PRINT  CUT  A  DEFICIT  INCURRED  BY  ANY  TOWN  ANC  THEN  ACCUMULATE  IT 
C 

60  WRITE  (OUT, 2000)  (NOCNAMI  INODE  +  LLL) ,LL=  1,31 ,DEF  ICT( IT0WN1 , I ITI ME 
TTLDF  (ITOWNI  *  TTLDF( ITOWN)  ♦  DEF  ICT( ITOWN) *TYMFAC 
OEFICTI  ITOWN)  =  0 
IFLAG  =  IFLAG  -  1 
IF  (IFLAG. EQ.O)  RETURN 
70  CONTINUE 
RETURN 


C 
C 


2000  FORMAT  (10X,3A4,«  HAD  A  DEFICIT  OP  •.FS^,*  MGD  IN  MONTH  NUMBER  • 
1,131 

END 
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SUBROUTINE  FILLST  ' 

C 
C 

C  REGULATES  THE  TRANSFER  OF  WATER  TO  AND  FROM  RESERVOIRS 

C 
C 

COMMON  IITIME,  ITIME,  MAXCAPC60),  NNCDE,  CI60),  STCAPYI60), 
1  T0WNI25),  TYMFAC 

CCMMCN  /M568/  IDVRZNI60),  IFILST(60)t  MAXMDCI25),  NODNAM ( 63, 3 ) , 
1  OUT,  SWSMAXI60.12) ,  TTLDFI25),  TTLDVI60),  TTLMDC125),  ICRDEPC60) 

CCMMCN  /S2356/  NPRINT,  STMAX (60, 12) ,  SWMIN(60,12) 

REAL  MAXCAP 

INTEGER  TOWN 

REAL  fAXMCC 

INTEGFR  CUT 
C 

c 

C  DIVERT  WATER  FROM  THE  FIVER  TO  A  RESERVOIR. 

C 

DC  20  JJ  =  l.NNODE 

IF  (ICRDER(JJI.EQ.O)  GO  TC  20 

INCDE  =  ICROERUJ» 

I  =  ICVRZN(INOCE) 

DIVERT  =  APIN1 I (MAXCAP(  INCDE)  -  STCAPY ( INCDE )) /TYMFAC, SWSMAXi INODE 
1, ITIME),  Q(I)  -  SWMIN( I, ITIME)) 

IF  (DIVERT. LE.O.)  GO  TC  20 

X  =  DIVERT  *  TYMFAC 

STCAPYC INODE)  =  STCAPYC INCOE )  ♦  X 
C 

C  WITH  PROPER  OPTICN  CCDE  PRINT  THE  AMOUNT  OF  WATER  DIVERTED 

C  FROM  THE  RIVER 

C 

IF  CCNPRINT/2)*2. EC. NPRINT)  WRITE(0UT,200C)  (NODNAM(  INODE, LL )  ,  LL 
1  =  1,2)  ♦  X  ,  IITIME 
C 

C         REMOVE  THE  DIVERSION  FRCM  THE  STREAMFLCW  AND  ACCUMULATE  THE 
C  DIVERSIONS 

C 

TTLDV(INCCE)  =  TTLDVC INCDE)  ♦  X 

DO  10  J  *  I, NNCDE 
10  QCJ)  =  Q(J)  -  DIVERT 
20  CONTINUE 
C 

C         TRANSFER  WATER  FROM  CNE  RESERVOIR  TO  ANOTHER 
C 

DO  30  INODE  =  i.NNODE 

IF  (IFILSTUNODE).EQ.O)  GC  TO  30 

I  =  IFILSTUNODE) 

XCHANG  =  AMIN1(STMAX(INCDE,ITIME)*TYMFAC, STCAPYC INODE ) , MAXCAPC I ) - 
lSTCAPY(I) ) 

STCAPY(INODE)  =  STCAPYC INCDE )  -  XCHANG 

STCAPY(I)  =  STCAPYC I)  ♦  XCHANG 
C 

C         WITH  THE  PROPER  OPTION  CODE  PRINT  THE  AMOUNT  OF  WATER  THAT  WAS 
C  TRANSFERRED  FRO*  CNE  RESERVOIR  TO  ANOTHER 

C 

IF  C(NPRINT/3)*3.EQ. NPRINT)  WRITE  COUT,2010)  XCHANG,  ( NODNAMC INODE 
1,LL),LL=1,3) ,(N0DNmM(I,LD,LL=1,3).IITIME 
30  CONTINUE 

RETURN 


C 

c 


2000  FORMAT  C81X,3A4,»  DIVERTEC  •,F6.1.»  MG  IN  MONTH  NUMBER  »,I3) 

2010  FORMAT  (5AX.F7.1,*  MG  WAS  DIVERTED  FROM  •,3A4,»  TO  f,3A4,«  IN  MONT 

1H  NUMBER  ■ ,13) 
END 
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SUBROLTINE    KEPTRK 

C 

c 

C  KEEPS  TRACK  OF  CERTAIN  STATISTICS  ON  A  MONTHLY  BASIS 

C 

c 

COMMON  IITIME,  ITIME,  MAXCAPI60),  NNCDE,  Q(60),  STCAPY(60). 
1  TCWNI25),  TYMFAC 

COMMON  /M78/  CLASS(7),  FREQ(10,8),  MAXST(60),  MINSTI60I,  NOBS, 
1  OBS(IO),  OBSMAX(IO),  OBSMIN(IO) 

CCMMCN  /R2/  ITRZFL,  ITRZST,  QTR ACE ( 1 20,9 ) ,  FLTRAZ ( 12C ,9)  ,  IWANTX 

REAL  ^AXCAP 

INTEGER  TOWN 

REAL  *AXST,  VINST 

INTEGER  CBS 


C 
C 


DC  20  IOeS  =  1,N0BS 

INODE  =  OBS( IOBS) 

IF  (  ITRZFL. NE. 01  QTRACE ( I  IT  IMF ,  IOBS »  =  Q(INODE) 

IF  (I1RZST.NE.0)  FLTRAZ ( I  IT IME,  IOBS >  =  STCAPY  (  INODE  I 

OBSMAXI  IOBS)  =AMAX1(CBSMAX( I OBS ) , 0 ( I  NODE ) ) 

0BSVIN(  ICRS)  =AMIN1(CBSHIN(  IOES),0( INODE  1  ) 

DC  10  I  =  1.7 

IF  (Q(INODE).LT.CLASS(I  )  )  GO  TO  20 
10  CONTINUE 

I  =  8 
20  FREGCIOBStll  *  FREO(IOBS.I)  ♦  1. 

IF  (IWANTX. NE. 01  CALL  WANTX1 

DO  30  INODE  =  l.NNOOE 

MAXSTUNCCEJ  =AMAX  1  (  MAXSTC  INODE  ),  STCAPY  (  I  NODE  I) 
30  MINST(INODE)  =AMI Nl ( P I NST ( INOCE ) , STC APY ( I  NODE  II 

RETURN 


END 
SUBROUTINE  APRINT 


C 

C 

C         PRINTS  ALL  OF  THE  SUMMARY  RESULTS 

C 

C 

COMMON  IITIME,  ITIME,  MAXCAPI60),  NNCCE,  C(60),  STCAPYC60), 
1  TOWN (25),  TYMFAC 

CCMMON  /M78/  CLASS(7),  FREC(10,8),  MAXST(60),  MINSTI60I,  NOBS, 
1  OBS(IO),  OBSMAX(IO),  OBSMIN(IO) 

CCMMON  /M568/  IDVRZNC60),  IFILST160I,  MAXMDC(25),  NODNAMI 63, 31 , 
1  OUT,  SWSMAX(60,12>,  TTLDFI25),  TTLDV<60),  TTLMDC<25»,  I0RDERC60) 
COMMON  /M35/  DEFICTI25I,  DEMAND! 25, 12  I ,  GWMAX( 25,12 1  ,  IFLAG, 

1  IPRELM,  IRTURN(25),  MDCCAP(25),  NTIHE,  NTOWN,  PR IORY( 25, 16) , 

2  RTURNI25.12I,  SWMAX( 25  ,  12  I  •  WELCAPI 25, 10 ,12 ) 

COMMON  /M12/  DA(10),  IMPMAX(12I,  IMPMINC12),  IMPNOD,  IMPQ,  IMPREC, 

1  IN,  NRCRDS,  NWELL,  PERCNT (60  ,12 ) ,  00(60),  QSUBAC60),  REC0RD(60), 

2  RESDA(60),  RES0OB(60),  SUBADA(60),  WELN0D(60),  WELPCT(60) 
COMMON  /R2/  ITRZFL,  ITRZST,  OTRACF ( 120, 9) ,  FLTRAZ! 120,9 1 ,  IWANTX 
DIMENSION  ITITLEI4I,  DUPMY(120),  N0DE(4) 

OIMENSION  IVAR(4) 

DATA  M,  N2,  N3  /  4HFLCW, 4HST0R,4H      / 
REAL  POCCAP,IMPMAX,IMPMIN,IMP0 
INTEGER  PRIORY,  RECORC,  RESOOB, WELNOD 
REAL  PAXCAP 
INTEGER  TCWN 
REAL  *AXST,  MINST 
INTEGER  OBS 
REAL  PAXMDC 
INTEGER  CUT 

INTEGER  FCRMT(10),FRMAT(4),FRMATT(4> 
C 
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C         IGNORE  COMPILER  WARNINGS  

C 

DATA    IVAR    /    16HA  BCD  / 

$ 
IFY0  38I    SIZE    WRN.*********************************************************** 
DATA    FRMAT    /16HC1H+  X,3A4)  / 

$ 
IEY03  8  I    SI  ZE    WRN. ************************************************ *********** 
DATA    FRMATT    /16H(1H+  X,F9.2)  / 

$ 
IEY0  38I    SIZE    WRN.*********************************************************** 
DATA    FORMT    /    40H.       1,    14,    27,040,    53,    66,    79,    92,105,118         / 

$ 
IE Y038  I    S I ZE    WRN.*********************************************************** 
C 
C 

C  PRINT  A  HYCROLOGICAL  SUMMARY  SHEET  FOR  EACH  RIVER  POINT  DESTRED 

C  WITH  A  FREQUENCY  TABLE 

C 

WRITE  (OUT, 2000)  (CLASS(II,I  =  1,7) 
DC  3C  I  =  i.NOBS 
TCTAL  =  0 
DC  10  J  =  1,8 
10  TOTAL  =  TCTAL  ♦  FREO(I.J) 

DC  20  J  =  1,8 
20  FREQMtJI  =  FREO(I,J)/TCTAL 

LL  =  CBS( I) 
30  WRITE  (OUT, 2010)  (  NODNAMI LL , LLL ) ,LLL=1,  3 )  ,OBSMAX(I)  ,OBSMIN(I)  , 
1 (FREO(I ♦ J) ,J=  It  81 
WRITE  (OUT,  2020) 
L  =  1 
C 

C  PRINT  THE  PRIMARY  SUPMARY  SHEET 

C 

DC  50  I  =  1,NNCDE 
IF  (TCWN(L).NE.I)  GO  TO  40 

WRITE  (OUT, 2030)  ( NODN  M*  (  I,LL  ),LL=  1,3)  ,MAXCAP(I),  MAXST(I)  ,MINS 
1T(I),  TTLDV(I),  TTLMDC(L),  MAXMDC(L),  TTLDF(L) 
L  =  L  ♦  I 
GC  TO  50 
40  WRITE  (OUT, 2030)  ( NCDNAM  I  ,LL )  ,LL=  1,3)  ,MAXCAP(I),  MAXST(I)  ,MINS 

ITU),  TTLDV(I) 
50  CONTINUE 
C 

C         WITH  PROPER  OPTION  CODE  PRINT  A  FLCW  TRACE  AT  ANY  RIVER  PCINT 
C  DESIRED 

C 

IF  (  ITRZFL.LE.l)  GO  TC  90 
WRITE  (OUT, 2040) 
DO   60  I  =  I.NOBS 

J  =  cesd) 

FPMATI2)  =  F0RMT(I*1) 
60  WRITE  (OUT, FRMAT)  ( NOON AMI  J , K ) ,K=1 ,3  I 

DC  70  I  =  1,NTIME 

WRITE  (OUT, 20501  J 

DO   7C  J  =  1,N0BS 

FPMATT(2)  =  FORMTU  +  1) 
70  WRITE  (0UT,FRMATT1  QTPACE(I,J) 
C 

C         WITH  PROPER  OPTION  CODE  PRINT  A  STORAGE  TRACE  FOR  ANY  RIVER 
C  PCINT  OESIRED 

C 


90  IF  (ITRZST.LE.il  GO  TO  125 
WRITE  (OUT, 2060) 
DO  100  I  =  1,N0BS 
J  =  CES( I) 

FRMATI2)  =  FORMTCI+1) 
100  WRITE  (OUT, FRMAT)  (NODNAMC J,K I ,K=1 ,3 ) 
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00    110    I    =    l.KTIME 
WRITE    (OUT, 2050)     I 

00  11C    J    =    l,NOBS 
FRMATT(2)    =    FCPMT(J»1) 

110  WRITE  (OUT.FRMATT)  FLTRAZU.J) 
C 

C         WITH  PROPER  OPTION  CODE  CALL  THE  SPECIAL  PLCT  ROUTINE  »WANTX« 
r 

125  IF  (  UANTX.NE.CI  CALL  WANTX 
C 

C  WITH  PPCPER  CPTICN  CCCE  PLCT  A  HYCRCGRAPH  OF  ANY  FLOW  TRACE 

C  OR  STORAGE  TRACE  DESIRED, 

C 

IF  (  ITRZFL.EC.O.AND. ITRZST.EO.OI  RETURN 
130  READ  (IN, 1000)  ( I TI TLE(  I )  ,NODE ( I ) , I =  1,4) 

DG  300  I  =  1,120 
300  DUMf'YU  )  =  0 
J  =  1 

DC  140  I  =  1,4 

IF  (  ITITLEI  II.EQ.Nl)  J  =  J*2 
IF  (  ITITLE(I).E0.N2)  J  =  J*3 
IF  (  ITITLE(  II.EQ.N3)  GC  TC  150 
140  CCNTIfvUE 

1  =  5 

150  NUM  =1-1 

IF  (MjM.EQ.O)  RETURN 

DO  160  I  =  l.NOBS 

IF  (NCDE(1).EC.CBS( I) )  Jl  =  I 

IF  (NCDE(2).EQ.0BS(I) )  J2  =  I 

IF  (NCDE(3).EQ.0BS(I) )  J3  =  I 

IF  (KCDE(4).EQ.0BS(I))  J4  =  I 
160  CONTINUE 

WRITE  (OUT, 20701 

DC   170  I  =  l.KUM 

L  =  NCDE( I) 
170  WRITE  (OUT, 2080)  IVAR(I),  ITITLE(I),  ( NODNAM( L ,LL ) ,LL=  1,3) 

WRITE  (GUT, 2090) 
C 

C  CALL  PLOT  ROUTINE  •PLOTER'  WITH  THE  DESIRED  PARAMETERS 

C 

IF  <( J/6)*6.EQ.J)  GO  TO  ( 130 ,260, 130, 270  I ,NUM 

IF  (( J/3)*3.EQ.J)  GO  TO  (220,230,240,250),  NUM 

GC  TO  (180,190,200,210),  NUM 
180  CALL  PL0TER(QTRACE(1, Jl) , DUMMY, DUMPY, DUMMY) 

GO  TO  130 
190  CALL  PL0TEP(CTRACE(1,J1),CTRACE(1,J2), DUMMY, DUMMY) 

GO  TO  130 
200  CALL  PLCTEP(CTPACE(1,J1),QTRACE(1,J2I,QTRACE(1,J3I, DUMMY) 

GG  TO  130 
210  CALL  PL0TEP(0TRACE(l,Jl),0TRACE(l,J21,GTRACE(l,J3),QTRACE(l,J4) ) 

GO  TO  130 
220  CALL  PLOTER  (FLTRAZ ( 1  ,J1 )  .DUMMY, DUMMY, CUMPY) 

GO  TO  130 
230  CALL  PLOTER  (FLTRAZ ( 1 , J  1 )  .FLTRAZ ( 1 , J  2) .DUMMY, DUMMY ) 

GO  TO  130 
240  CALL  PLOTER  ( FLTRAZ ( 1 , J  1 ) , FLTRAZ( 1, J  2 ) , FLTRAZ( 1 , J3 ) .OUMMV) 

GO  TO  130 
250  CALL  PLOTER  (FLTRAZ ( 1  ,J  1 ), FLTRAZ ( 1 ,J2 ), FLTRAZ( 1 ,J3 ) ,FLTRAZ( I , J4 ) ) 

GO  TO  130 
260  CALL  PL0TEP(CTRACE(1,J1),FLTRAZ(1,J2), DUMPY, DUMMY) 

GC  TC  130 
270  CALL  PL0TER(QTRACE(1.J1),FLTRAZ(1,J2).QTRACE(1,J3),FLTRAZ(1,J4) ) 

GC  TO  130 
C 

c 

1000  FCRMAT  (4(A4,I2)) 

2000   FCRMAT  (1H1///  4X,  'CBSEPVATION  NODE* ,3X, 'MAX  MONTHLY  FL0W»,3X,»M 
UN  MONTHLY  FLOW* ,25X, • FREGUENCY  TABLE   (MGD)»  /  28X, • ( MGD) • , 14X, • ( 
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2MGD)  i,7X«0.,»7C    <0<     *,F4.0),     •    <0 «     I  ~—     

2010    FORMAT    (6X,3A4,7X,F8.0,13X,F6.2,6X,8(4X,F5.3)  ) 

2020  FORMAT  (1H1///  2X, 'RIVER  PCI  NT* ,3X, • STORAGE  CAPACITY • ,2 X  ,* MAX  STOR 
1AGF* ,2X,*MIN  STORAGE*, 2X, *FLOW  DIVERTED* , 2X, • WATER  PURCHASED •  ,2X , » 
2MAX  MONTHLY  PURCHASE*  ,2X , 'TOTAL  DEFICIT*  /  22X, • ( MG )  • , 12X, • ( MG  )  • 
3,9X,*(MG) ',9X, MMG) ',13X, *(MG)* , 15X,*(MG) ', 14 X,* (MG)  •   » 

2030  FORMAT  ( 2X ,3A4 , 6X,  F6 . 0, 10 X, F6 .0, 7X, F6.0t 5X,F 8.0 , 1 IX, F 7. 1 , 13X ,F6.  1, 
110XfF7.0l 

2040  FCRMAT  ( 1H I// ,40X, • FLOW  TRACE  AT  INDICATED  PCINT    (MGD)*  //5Xt 

1  « MONTH*  I 

2050  FORMAT  l/,6X,I3) 

2060  FORMAT  ( 1H1//, 30X, 'RESERVOIR  STORAGE  TRACE  AT  INDICATED  POINT    (M 

1G  > •   //5X,  *MONTH*  ) 
2070  FORMAT  (1H1  ////  10X,*A  PLCT  FOR  THE  PERIOD  OF  RECORD*  //20X,'IF  M 

10RE  THAN  ONE  SCALE  APPEARS,  THE  FIRST  IS  FOR  SYMBOLS  A  ANO  C*  / 

2  20X,  *AND  THE  INNER  IS  FCR  SYMBCLS  8  AND  C.   »///30X*X  REPRESENTS 

3  MORE  THAN  ONE  VARIABLE  AT  A  POI.NT*   I 

2080  FORMAT  (30X,A1,»  REPRESENTS  • , A4, *  FOR  •  ,3A4) 
2090  FCRMAT  (////) 
END 

SUBRCLTINE  PLOTER  (A,B,C,C) 

C 

c 

C  PLCTS  A  TIME  SERIES 

C 

C 

COMMON  /M568/  IDVRZNC60I,  IFILST(60>,  MAXMDC(25),  NOCNAM (63, 3 ) , 
1  OUT 

CCMMCN  /K35/  DEFICTI25),  CENAND( 25, 12 ) ,  GWMAX ( 2 5 , 12 > ,  IFLAG, 
1  IPRELM,  IRTURNI25),  MDCCAP(25),  NTIME 

DIMENSION  A(120),  B(120),  C(120),  D(120),  DUMMYK1),  X(13),  DUMMY2 
1(1),  Y<13),  DELTAI2I,  RANGE(2),  SCALE(2,13>,  IDUMMY(l),  LINE(125) 
DIMENSION  MONTH! 11) 

DATA  KONTH  /4HJAN  ,4HFEe  ,4HMAR  ,4HAPR  ,4HMAY  ,4HJUNE ,4H JULY, 4HAUG 
1  ,4HSEPT,4H0CT  ,4HN0V   / 

DATA  CUMMY1,  X,  0UMMY2,  Y  /  40000. ,20000 ., 18750. , 12500. , 6250. , 
13125., 2500., 1875., 12 50.,  6 25., 500., 250., 12 5., 0., 400., 200., 150., 100 „ 
2,  50. ,25. ,20. ,15. ,10. , 5. ,4. ,2. ,1. ,0.   / 
DATA  IA, IB,IC,  ID,  IX, II,  IBLANK  /1HA , 1HB, 1HC , IHD, 1HX, 1HI , 1H  / 
INTEGER  CUT 
C 
C 

C  FIND  MAXIMUM  VALUE  OF  EACH  SERIES 

C 

AAMAX  *  0 
BMAX  =  0 
CMAX  =  0 
DNAX  =  0 

DC  10  I  =  1, NTIME 
IF  (AIU.GT.  AAMAX1  AAMAX  =  A(I) 
IF  (B(I).GT.  BMAX)  BMAX  =  B(I) 
IF  (C(I).GT.  CMAX)  CMAX  =  C t 1 1 
IF  (0(1). GT.  DMAX)  DMAX  =  D<I) 
10  CONTINUE 
C 

C         COMPUTE  SCALE  FACTORS 
C 


RANGE(i)  =  IFIX(AMAX1 (AAMAX, CMAX)  ♦  .999) 

RANGE(2)  =  IFIX(AMAX1(  BMAX, DMAX)  ♦  .999) 

DC  30  I  =  1,2 

DO  20  J  =  1,12 

IF  (RANGE(  D.GT.Xt  J)  )  GO  TO  30 
20  CONTINUE 

J  =  13 
30  DELTA(I)  =  Y(J-l) 
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40 


50 

60 


K  =  0 

IF  (0ELTA( 1) .LE.Y(J-2).ANC.DELTA(1) .GE.Y(JI) 

IF  (K.EQ.O)  GO  TO  40 

DELTA  (I)  =AMAX1(CELTA(1 ),DELTA(2I I 

DELTA(2)  *  DELTA(l) 

DO  60  I  =  1,2 

SCALF (1,11  =  0 

SCALE(I,13»  *  DELTAC  I)*125. 

DC  50  J  =  2,  12 

SCALE(I,J)  =  DELTA(I)*10.*FLOAT(J-1I 

WRITE  (OUT, 2000)  ( SCALE ( I  ,  J) , J=i , 13) 

WRITE  (CUT.2010) 

WRITF  (OUT, 2040) 
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C 
C 
C 


70 


80 
SO 

100 

110 

120 
130 


140 
150 


C 
C 


K  =  1 


PLCT  EACH  OF  FOUR  POSSIBLE  TIME  SERIES  EACH  MONTH 


DO  15 

12  = 

13  = 

14  = 

J  =  0 

LINE( 

DO  70 

LINE( 

II  = 

IF  (B 

IF  (C 

IF  (D 

LINE( 

IF  (L 

LINE( 

IF  (L 

LINE( 

IF  (L 

LINEC 

GO  TO 

LINEi 

GO  TO 

LINE! 

GO  TO 

LINEC 

IF  CL 

IF  CL 

WRITE 

L  =  I 

IF  <L 

WRITE 

GO  TO 

WRITE 

CCNTI 

WRITE 

WRITE 

WRITE 

WRITE 

RETUR 


=  l.NTIPE 


C  I 
C 
C 
C 


J)  =  IBLANK 
J  =  1,125 
J)  =  IBLANK 
MII/DELTM1I  +  1.5 
MAX.NE.O)  12  =  B(I  )/DELTA(2) 
MX.NE.01  13  =  C(  I  )/DELTA(l) 
fAX.NE.O)  14  =  C(I >/DELTA(2) 
III  =  IA 

INE(I2).NE.IBLANKI 
121  =  IB 
INF(  m.NE.IBLANKI 

13)  ■  IC 
INE(  14). NE. IBLANK) 

14)  =  IC 
130 

IX 


♦  1.5 

♦  1.5 

+  1.5 


GC  TO  100 
GO  TO  110 


GC  TC  120 


121  = 
80 

13)  »  IX 
90 

14)  =  IX 

INECD.EQ.  IBLANK)  LINE(l)  =  II 
INE(125).E0. IBLANK)  LINEC125)  =  II 

(OUT, 20501 

-  (1/121*12 
.EO.O)  GO  TO  140 

(OUT, 2020)  M0NTMl).(LINE(J),J=l,125l 

150 

(OUT, 2030) 
NUE 

(OUT, 2040) 

(CUT, 2010) 

(OUT, 2000) 

(OUT.2000I 
N 


I.  (L  INEC J),J=1,125) 


(SCALE(2,J),J=1,13) 
(SCALE(1,J) ,J=1,13) 


2000 
2010 
2020 
2030 
2040 
2050 


FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
END 


(      3X,  12(F6.0,4X),3X,F6.0) 
(7X,12(»*« ,9X) ,4X,»*» I 
(4X.A3,125A1) 
(4X,I3,125A1) 
(7X,1HI,123(1H*),1HI) 
(7X,1HI,123X,1HI ) 
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C 
C 

C  PRINTS  A  LIST  OF  THE  INPUT  DATA 

C 
C 

COMMON  IITIME,  ITIME,  MAXCAPI60),  NNOOE,  C(60),  STCAPY(60), 
1  TCWN(25),  TYMFAC 

COMMON  /t*35/    DEFICT(25),  CEMANO(  25,  12  )  ,  GWMAX (25, 12 ) ,  IFLAG, 

1  IPRELM,  IRTURNI25),  MDCCAP(25),  NTIME,  NTOWN,  PR IORY( 25  ,16 ) , 

2  RTURN<25,12),  SWMAX ( 25 , 12 ) ,  WELCAP( 25, 10, 12 ) 

COMMON  /M568/  IDVRZNI60I,  IFILSTI60),  MAXMDC(25),  NODNAM ( 63, 3 ) , 
1  OUT,  SWSMAX(60,12),  TTLDF<25),  TTLDV(60),  TTLMDC(25I,  ICRDER(60) 

COMMCN  /M24/  EFF(60,12),  UNIT(60,12),  WELL(25,10) 

CCMMCN  /S2356/  NPRINT,  STPAX ( 60 ,1 2 ) ,  SWM I N (60 , 12 ) 

RE*L  MCCCAP 

INTEGER  PRIORY 

INTEGER  CUT 

INTEGER  WELL 

INTFGER  TCWN,  TITLE 
C 
C 

WRITE  (OUT, 2050) 

IRES  =  0 

DC  RO  I  =  1, NTOWN 

J  =  TCWN(  I) 

WRITE  (OUT, 2000)  (  NOCls  AM  J,K  )  ,K  =  1  ,3  ) 

KCUNT  =  1 

GC  TC  90 
C 

C         PRINT  THE  SOURCES  FOP  EACH  TOWN 
C 

5  DO  57  J  =  1,8 

L  =  PRIORYU  ,J) 

IF  (L.EO.O)  GO  TO  60 

IF  (L-2)  10,30,40 
C 

C  PRINT  THE  CAPACITY  OF  A  WELL 

C 

10  DO  20  M  =  1,10 

IF  (WELLd.MI.EQ.OI  GC  TO  57 

WRITE  (OUT, 2010)  WELL ( I ,M  )  ,WELCAP( I , M, 1 ) 

KOUNT  =  2 

GO  TO  90 
15  CONTINUE 
20  CONTINUE 

GC  TO  57 
C 

C  PRINT  THE  CAPACITY  OF  A  STREAM  DIVERSION 

C 

30  WRITE  (OUT, 20201  SWMAXM.l) 

GO  TO  50 
C 

C  PRINT  THE  YIELD  CAPACITY  OF  A  RESERVOIR 

C 

40  K  *  L  -  10 

WRITE  (OUT, 20301  CNOONAMf K,KK ) ,KK  =  1,3) , STMAX( K , 1) 
50  KCUNT  =  3 

GO  TO  90 
53  CONTINUE 
57  CONTINUE 
60  J  =  J  +  1 

IF  (J.GT.16)  GO  TO  80 

L  ■  PRIORYU, J) 

IF  (L.EO.O)  GO  TO  60 

IF  (L  -  2)  70,10,40 
C 
C  PRINT  THE  CAPACITY  OF  A  CONNECTION  TO  AN  EXTERNAL  SOURCE 

40 
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70    WRITE     (0UT,204C)     MCCCAP(I) 

KCUNT  =  4 
90  IPES  =  IRES   ♦  1 

IF  (IRES.GT.NNCDE)  GC  TC  (5, 1 5 ,53,60  ) ,K0UNT 
C 

C  PRINT  EACH  RESERVOIR  INVOLVED  IN  TRANSFER  OF  WATER  ANC  THE 

C  CAPACITY  OF  THE  PUMP  OR  PIPELINE  THAT  TRANSFERS  THE  WATER 

C 

IC  =  IDVRZNi  IRES! 
I FIL  =  IFILST(IRES) 
IF  f  I  CEO. 01  GO  TO  120 
IF  (IFIL.EC.O)  GC  TO  100 

WRITE  (CUT, 2070)  ( NCDNAM IRES,KK) ,KK=1 ,3 ) , (NCDNAM ( I D,KK) ,KK= 1, 3  ) , 
1  SWSMAXt IRES.l I, (N0DNAM(IFIL,KK),KK=1,3) , STMA X( I  RES, I) 
GC  TO  130 
100  WRITE  (OUT, 2080)  ( NODNAM ( IRE S ,KK) ,KK=1 ,3 ) , ( NCDNAM ( ID, KK ) , KK= 1 , 3  )  , 
1  SWSMAX( IRES,1) 
GC  TO  130 
120  IF  ( IFIL.EQ.O)  GO  TO  90 

WRITE  (OUT, 20^0)  (NOON AM IRES,KK ) ,KK= 1 , 3 ) , ( NODNAM ( IF IL , KK ) ,KK=1 , 3 ) 
1  ,  STMAX( IRES,1) 
130  GO  TO  (5,15.53,60),KCINT 
80  CCNTINUE 

WRITE  COUT, 20601 
RETURN 
C 
C 
2000  FORMAT  (  IX, 3A4) 

2010  FCRMAT  (16X,'WELL  NUMeEP  ■  ,  I  2, 3X, F6. 2 ) 
2020  FORMAT  (16X, 'STREAM  DIVERSION  »,F6.2) 
2030  FORMAT  ( 16X, 3A4, 5X, F6  .2  ) 
2040  FCRMAT  < 16X, • TREATED  IMPORT  «,2X,F6.2» 

2050  FORMAT  (/////  AX, 'LIST  CF  TOWNS  WATER  RESOURCE  SUPPL  IES •  ,  15X, 'L  1ST 
1  OF  RESERVOIR  INFLOWS  AND  OUTFLOWS'///  IX, 'TOWNS ', 10X, • SOURCES' , 8X 
2, 'PUMP  CAPACITY'  /  35X,MMGD)»   I 
2060  FCRMAT  ( 1H1) 
2070  FCRMAT  ( 1H*. 46X ,3A4, '  DIVERTS  AT  ',3AA,'  AT  »,F5.1,'  MGD  AND  FILLS 

I  ',3A4,'  AT  ',F5.1t«  MGC  '1 
2080  FCRMAT  ( 1H+.46X  ,3A4, •  DIVERTS  AT  ',3A4,«  AT  ',F5.1,'  MGD') 
2090  FORMAT  (  1H+, 46X ,3A4, •  FILLS  ',3AA,'  AT  ',F5.1,'  MGD') 
END 

SUBPCLTINE  WANTX 
C 

C         SETS  UP  PARAMETERS  FCR  A  SPECIAL  PLCT  TO  BE  CCCED  AND  THEN  CALLS 
C  THE  PLOT  ROUTINE  'PLOTER' 

C 
C 

CCMMON  IITIME,  ITIME,  MAXCAP(fcO),  NNODE,  0(60),  STCAPY(60), 
1  T0WN(25),  TYMFAC 

CCMMON  /M12/  DA(IO),  IMFMAX(12),  IMPMIN(12),  IMPNOD,  IMPC,  IMPREC, 

1  IN,  NRCRDS,  NWELL,  PERCNT( 60 ,12 ) ,  00(60),  0SUBA(60),  REC0RD(60), 

2  RESDA(60),  RES008(60),  SUBA0A(60),  WELNOD(60),  WELPCT(60) 
COMMON  /M35/  0EFICT(25),  DEMAND( 25 , 12 ) ,  GWMAX (25, 12) ,  IFLAG, 

1  IPRELM,  IRTURN(25),  MDCCAP(25),  NTIME,  NTOWN,  PRIORY( 25 , 16) , 

2  RTURN(25,12),  SWMAX ( 25  ,  12 ) ,  WELCAP (25, 10, 12 ) 

COMMON  /M78/  CLASS(7),  FREQ(10,8),  MAXST(60),  MINST(60),  NOBS, 
1  OBS(IO),  OBSMAX(IO),  CBSMIN(IO) 

COMMON  /M568/  IDVRZN(60),  IFILST(60),  MAXMDC(25),  NODNAM (63, 3 ), 
1  OUT,  SWSMAX(60,12).  TTLDF(25),  TTLDV(60),  TTLMDC(25),  ICRDER(60) 

CCMMCN  /M24/  EFF(60,12),  UNIT(60,12),  WELL(25,10) 

CCMMON  /S2356/  NPRINT,  STVAX (60,12 ) ,  SWMIN(60,12) 

COMMON  /R2/  ITRZFL,  ITRZST,  QTRACE ( 120,9) ,  FLTRAZ( 120,9) ,  IWANTX 

REAL  IMPMAX,  IMPMIN,  IMPQ 

REAL  MDCCAP 

REAL  MAXST,  MINST  • 
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REAL  MXMDC 
REAL  VAXCAP 

INTEGER  RECORD,  RESOOB*  WELNOD 
INTEGER  PRIORY 
INTEGER  OBS 
INTEGER  OUT 
INTEGER  TCWN 
INTEGER  WELL 

DIMENSION   AI120),  B(12C),  C<120),  0(120) 
C 

c 

C         CALL  PLCT  ROUTINE  'PLOTER* 
C 

WRITE  (OUT, 2000) 

CALL  PLOTEP  (A,B,C,C) 

RETURN 

ENTRY  WANTX1 

IF  (  IWANTX.EQ. 98651)  GO  TO  1C 

IWANTX  =  98651 

DO  5  I  =  1,120 

A(I)  =  0. 

B(I)  =  0. 

C(I)  =  0. 
5  D(I)  =  0. 
C 

C  USE  THE  FOUR  TIME  SERIES  •A«,«B,,,C,t  AND  'D*  AND  CODE  ANY 

C  VARIABLE  DESIRED  TO  BE  PLOTTED 

C 

10  A(IITIME)  =  0. 

C(IITIME)  -    0. 

D(IITIME)  =  0. 

B(IITIME)  =  0. 

RETURN 

2000  FORMAT  (  1H1////5 X, • AN  EXTRA  PLOT  FOR  THAT  SOMETHING  SPECIAL  YOU  WA 
1NTED«  /////> 
END 


C 
C 
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SAMPLE  INPUT  DECK 


The  following  is  a  list  of  the  data  cards  which  causes  the  results  of  Appendix  IV. 
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COLUMN 

12  3  4  5  6  7  0 

1  734  5678^0 123456 7 6?0123 4 567 8 9012345678901234567890123 4 5678901234567890123456769U 

^6 

EXAMPLE  RON  FOR  REPOR *  -  ONE  FEASIBLE  ALTERNATIVE  FOR  1995 

48117  20   9  38   3  24210   2   2   1 
BURLINGTON 
V4PLE  MEADOW 
WILMINGTON 
READING  WELL 
READING 
N  READING 
BEAR  MEADOW 
4ARTIN  BROCK 
EISENHOUR  PD 
LYNN FIELD 
L  Y  N  M 

^AUGUS  RTVER 
H  A. W  K  r  S  POND 
W  A  L  r' r  M  POND 
RI RCH-BREE^S 
PFABODY 
SUNTAUG  LAKF 
S  p  R  T M  G  LAKE 
S  MIDDLETON 
PEABODY  DIV 
MORRIS  BR  00^ 
DA.NVERS-MIDD 
DAM VERS  WELL 
MIDDLETON  PD 
EMERSON  POND 
BOSTON  BROOK 
MICHOL  BROOK 
BOXFORD 
SALFM-BEVRLY 
WFNHAM 
HAMILTON 
TOPSFIELD 
GEORGETOWN 
I PSWI CH  GAGF 
IPSWICH  WELL 
ROWLEY 

IPSWICH  WELL 
R  c s  =  30-B 
WEMHAM  LAKE 
IPSWICH 
ESSEX 

GLOUCESTER 
WALLACE— DIKE 
R ARSON-GOOSE 
KLONDIKE 
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COLUMN 

12  3  4  5  6  7  d 

123 4 567890123456789 0123 4 567 8 9 01234567890123456 789012345 6 78901 2 34567o90l23456 7 d9u 

MANCHESTER 

ROCKPORT 

IPSWICH  DAM 

IMP0RT01  3   3 

RcC09noi  11   111111111111111112   2^22 

^r^Q^02  2222222   2   22222222222   2   2222 

WFLNOP01  2   2   2   3   3   4   5   6   8   8   8   8   d  1 0  1 U  2 1  Z 1  ZZ     <L 3  2d  30  3  u  3o  js  i 

WELNOD02  32  33  33  35  36  37  38  38  39  39  41  41  46  47 

TOWN   01  1   3   5   6  10  11  16  22  28  29  3u  31  32  33  j6  4u  41  4^  46  47 

0DS    01  3   6  1)  28  29  34  37  38  48 

IRTURN01  48   3   5   6  10  4d  16  22  46  48  30  31  3  <i  48  46  40  4o  4o  4o  4o 

PERCNT01  067  067  067  067  067  067  067  067  Ob7  067  06  7  067 

PFRCNT02  132  132  13?  132  132  132  132  132  132  132  132  132 

PERCNT0  3  130  130  130  130  130  130  13u  130  130  130  130  13u 

DFRCNT04  012  0  12  012  012  012  012  012  012  012  012  012  012 

PrRCNT0  5  002  002  0  02  002  00  2  0  0 2  00^    0  0 2  002  002  00  Z    uO^ 

^"rRCNT06  02^  023  02^     023  02  3  023  02  3  023  023  023  023  023  02  3 

nrn^N1T07  127  i?7  i2.-f     127  127  127  127  127  127  12  7  127  127 

ppRCNTOR  3  37  337  337  337  337  3  37  337  337  3  37  337  337  3  37 

PFRCNT09  n?7  027  0270)27  027  027  027  027  027  027  027  027 

PERCNT10  06?  062  CS2  062  06?  06?  062  062  062  062  06Z    062 

PFRCNT11 

PCRCNT12 

PERCNT13 

PERCNT14 

PFRCNT15 

PFRCNT16 

DFRCNT17 

PFRCNT18 

P-RCNT19  035  035  035  0  35  035  035  035  0^5  035  0^5  0  35  0^3 

PFRCNT20  001  001  001  001  001  001  001  001  001  001  001  001 

PERCNT21  082  082  0  82  Od2  082  082  Od^  08  2  06  Z     ObZ    08^:  032 

PERCNT22  005  005  005  005-  005  005  005  005  005  005  005  005 

PERCNT23  006  006  006  006  006  006  006  006  006  006  006  006 

nrRrN'T?4  o21  021  021  021  021  021  021  021  021  021  021  021 

PERCMT25  OAO  040  040  040  040  040  040  040  040  040  040  040 

PFRCMT26  219  219  21Q  2 1 c'  219  219  219  219  219  219  219  219 

PERCNT27  054  054  054  054  034  054  054  034  034  054  034  034 

PFRCMT28  233  233  2  3  ">  233  233  23^  233  233  233  233  233  233 

PERCNT29  0  500050 00 50005 000500 0  50005 000500050005 00050UUDO 

DCRCNT30  034  034  034  034  034  034  034  034  034  034  0^4  0^4 

PFPCNT31  046  046  046  046  046  046  046  046  046  046  046  046 

PFRCNT32  175  175  175  175  175  175  175  175  175  175  175  175 

DERCNT33  023  023  023  023  023  023  023  023  023  023  023  023 

PFRCNT34  01 00 01000 10 0  01000100 010 0010001000100 01000100010 

PFRCNT35  005  005  005  005  005  005  005  005  005  005  005  005 

PFRCNT36  03  0  00300030003  0  0  0300030  0  0  300030  00300030  0  0300030 

PERCNT37  005  005  005  005  005  005  005  005  005  005  005  005 
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PERCNT38  005  005  005  005  005  005  005  005  005  005  005  005 

PFRCNT39  -225  225  225  225  225  225  225  225  225  225  225  ^25 

PERCNT40  034  034  034  034  034  034  034  034  034  0^4  034  034 

PERCNT41  004  004  004  004  004  U04  004  004  004  004  004  004 

PERCNT42 

P,rRCMT43 

PERCNT44 

PERCNT45 

PERCNT46 

PERCNT47 

orRCNT48 

RE5DA  01  io4 

RFSDA  02      1913     17  87  5  lb 

RESH.A  03     2^  9 

RFSDA  04  1^     59       9  61035       1       1       o 

SUPADA01     29     b~>  56       5       1      10     55    146     31     27 

SURADA02  15       1     bo      4       o  17 

SUCADA03     32    175     43    106     40     27      37    140     Id       8       4     24 

SUBADA04      4      4    180     27      3 

TFILST01  13  14  15         lb  ^4 

IFILST02  24  39 

MAXCAP01  .  2ubu 

MAXCAP02   2800  20000  19500  6000   5000  1500  4u0b 

MAXGAP03   3  0  0  0  17  0  0  0 

MAXCAP04        100000  12550    600  4000   3500   5000   il50   1360   ioOO 

5WM INFO  102 

SWMAX  01 

SWMAX  0  2 

SWMAX  03 

SWMAX  04 

SWMAX  OR 

SWMAX  06 

SWMAX  07  .   .. 

SWMAX  08 

SWMAX  09 

SWMAX  10 

SWMAX  11 

SWMAX  1? 

SWMAX' 13 

SWMAX  14 

SWMAX  15 

SWMAX  16 

SWMAX  17 

SWMAX  18 

SWMAX  10 

.SWMAX  2  0 

SWS^AXOl 

SWSMAX02 
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SWSMAX03 

SWSMAX04  • 

SWS^AX()5 

5WSMAX06 

SWSMAX07 

SWSMAX0  8 

SWSMAX09 

SWSMAX10 

sw^yax 1 l 

5W.SVAX12 

SWS^AXl-3 

TOO 

1  0  0 

IOC 

100 

100 

S  hi?  ^  AX  14 

100 

100 

1  00 

100 

100 

SWSVAX15 

SWSMAX16 

SWSMAX17 

200 

200 

200 

200 

200 

SWSMAX18 

S  \\'  S  K  A  X  1  9 

c;  i.l  r.  \>  A  X  P  D 

S'..'S"AX?1 

5WSMAX2  2 

SWSMAX2  3 

SWSMAX24 

SW5MAX2  5 

S'J.SMAX.26 

SWSVAX27 

SWS^AX28 

SWSMAX29 

250 

2  50 

250 

250 

250 

SWSMAX30 

SWSMAX31 

SWSMAX3  2 

SWSMAX33 

SWS.MAX3  4 

SWSMAX.35 

SWSMAX36 

SW5MAX3  7 

SWSMAX38 

400 

400 

400 

400 

40  0 

SWSMAX39 

250 

250 

250 

250 

250 

5WSMAX40 

SWSMAX41 

SWSMAX42 

SWSMAX43 

SW5MAX44 

SWS^AX4  5 

SWSMAX46 

S.WSMAX47 

SWSMAX48 

STMAX  01 

STMAX  02 

luu 


6.  ul; 


25o 


4Uu  4uU 

25J 
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ST  MAX 

03 

STMAX 

04  . 

S  T  M  A  X 

0  5 

S  T  M  A  X 

06 

STMAX 

07 

5  T  M  A  X 

08 

S  T  M  A  X 

09 

S  T  M  A  X 

10 

ST  MAX 

11 

ST  MAX 

12 

60 

60 

60 

60 

60 

6U 

6(J 

6  0 

60 

6o 

DO 

6u 

S  T  M  A  X 

13 

200 

200 

200 

200 

200 

200 

200 

ZUJ 

2ou 

10^ 

dUU 

<£O0 

STMAX 

14 

200 

200 

200 

200 

200 

200 

2  00 

200 

2<ju 

IvJ 

^OU 

<^  u  u 

ST  MAX 

15 

240 

240 

24U 

240 

^40 

240 

240 

Z40 

240 

24 

L^tU 

2  4u 

S  T  M  A  X 

16 

ST'' AX 

17 

^6 

5  6 

5  6 

5  6 

5  6 

56 

56 

56 

5o 

56 

56 

56 

c  T  M  A  X 

10 

1  ^1 

i  n 

1  2  1 

121 

121 

121 

121 

121 

121 

121 

121 

121 

S  T'MX 

19 

ST  MAX 

2  0 

STMAX 

2] 

ST  MAX 

22 

p 

2 

2 

2 

2 

2 

2 

2 

2 

z 

<L 

2 

ST  MAX 

23 

S  T  M  A  X 

2  4 

43 

43 

43 

43 

43 

43 

43 

43 

4J 

43 

^3  . 

43 

S  T  M  A  X 

25 

1  0 

10 

10 

10 

10 

10 

10 

10 

10 

10 

lo 

10 

STMAX 

26 

STMAX 

27 

STMAX 

28 

STMAX 

29 

250 

250 

25o 

2  50 

^50 

2^u 

25o 

ibo 

C'J\J 

25o 

t.  5  y 

Zbw 

STMAX 

3  0 

STMAX 

31 

STMAX 

32 

STMAX 

33 

STMAX 

34 

STMAX 

35 

S  T  M  A  X 

36 

S  T  M  A  X 

37 

STMAX 

38 

100 

100 

100 

100 

100 

100 

100 

loo 

100 

100 

iOO 

x  oo 

STMAX 

30 

240 

240 

240 

240 

240 

240 

240 

^4j 

240 

240 

^4u 

240 

STMAX 

40 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

S  T  M  A  X 

,41 

STMAX 

42 

50 

50 

50 

50 

50 

50 

50 

50 

5u 

5o 

5u 

5u 

S  T  M  A  X 

4  3 

45 

^+5 

45 

4  5 

45 

45 

45 

45 

45 

45 

45 

45 

STMAX 

44 

4  5 

45 

45 

45 

45 

45 

45 

45 

45 

4^ 

45 

42 

STMAX 

45 

8 

8 

8 

8 

8 

8 

8 

8 

a 
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8 

8 

STMAX 

4  6 

14 

14 

14 

14 

14 

14 

14 

14 

14 

14 

14 

14 

STMAX 

47 

10 

10 

10 

10 

10 

10 

10 

10 

10 

10 

10 

10 

S'T  MAX 

48 

RPSOOB01 

1 

1   1 

1 

1   1 

1 

RPT-OOR-02 

1 

1 

1 

1   1 

1    1 

1 
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D  E  M  A  N 

f)FMAM 
D  F  M  A  M 

[)  P  M  A  (\i 

0  E  M  A  N 
D  E  M  A  N 
D  E  M  A  N 

DEMAN 

r>  p  m  ^  po 

n  c  M  £  \J 
DFMAN 
DEMAN 
TEMAN 
■D  E  M  A  N 
DEMAN 
OF  MAN 
DEMAN 
n  F  m  A  N 
pi  P  y  A  M 
D  F  M  A  N 
RTURN 
RTURN 
R  T  U  R  N 

rturn 
R  T  U  R  N 
RTURN 
RTURN 
RTURN 
RTURN 
RTURN 
RTURN 
RTURN 
R  T  U  R  N 
RTURN 
RTURN 
R  T  U  R  M 
RTURN 
RTURN 
R  T  U  R  N 
RTURN 
MDCCA 
MDCCA 
PR  I  OR 
PRIOR 
PR  IOR 
PRIOR 
PRIOR 
PRIOR 


D01 
P02  ■ 
H>03 
004 
D05 
D06 
DO  7 
DO  8 

Din 
mi 

D12 
D13 
D14 
D15 
016 
D17 

ni  8 
D1  9 
D?0 

01 

0  2 

03 

04 

0  5 

06 

07 

08 

09 

10 

111 

121 

131 

14 

15 

16 

17 

18 
■  IP 

20 
P01 
P02 
Y01 
Y02 
Y03 
Y04 
Y05 
Y06 


771 
835 
78  } 
74  > 
799 
833 
7  69 
7^7 
7  84 
6  79 
538 
70  7 
734 
777 
720 
5  39 
779 
64Q 
615 

8 
2 


672 
730 
883 
635 
727 
7  16 
74  b 
659 
76  6 
62^ 
724 
659 
724 
729 
682 

4  6  4 
69  2 
595 

5  50 

3 
2 


?  5  3 
8  31. 
920 
721 
80  7 
753 
836 
710 
7  76 
69  3 
653 

710 
752 
82  3 
80^ 
50? 
721 
5  96 
578 


788 
764 
^69 

'46 
002 
843 
767 
718 
78  5 
755 
75  1 
718 
702 
8  72 
860 

6  90 

7  64 
66  3 


892 

854 

826 

9461 

843 

8  56 

84  5 

92  5 

812 

871 

7  80 
92  5 
812 

8  17 
8341 
!■)  3  1 
80  8 

8  1  8  1 
8281 

64 
16 


9  3  7 

9  54 
94  7 
012 
890 
9  6  0 
912 
9  77 
891 
8  0  8 
93  8 
977 
89  8 
982 
007 
2  10 
917 
00  5 
180 


16 


08 
U3 


944 

974 

1090 

1092 

94  5 

9  56 

94  2 

115  7 

9  4  8 

1173 

1  227 

1  157 

96  3 

95b 

103  5 

1 4  6  3 

98  6 

1277 

133  a 

16 
04 


993 

934 

1102 

906 

965 

9  26 

977 

10  8  6 

971 

120  7 

1149 

10ti6 

929 

9  59 

1083 

120  7 

10  3  7 

]  266 

12  96 

16 
0  4 


882 

819 
851 
897 
86  2 
8j8 
640 
881 
872 

104R 
8  74 
881 
871 
804 
843 

1051 
8  95 

10  04 
966 


08 


847 
76  1 
62  5 
768 
8  3  7 
o  u  3 
6  j4 
736 
b  16 

6  18 
762 

7  66 
6  ci4 
713 
74  2 
66  6 
6  28 
776 
742 

56 

14 


762 

741 
626 
726 
13j 
7j4 
7j2 
724 
790 
6  59 
761 
724 
8  70 
796 
64j 
6^:4 
776 
666 


759 
779 
7  64 
7  77 
7  78 
771 
779 
712 
782 
660 
76^ 
712 
o  52 
764 
711 
636 
791 
66  V 
6  4  9 


04   04   08   14 


72   8 
16   2 

16   2 


3  5   35   35   35   28   14   07   07   14   245  315  35 


8 
8 


u 

4 
4 


2 
2 

2 


2 
2 
2 


4 
4 
4 


7  9  1 
7  9  1 
7    9   1 


35   45 


40   40 


1  48 

1  48 

1  48 

1  48 

25  48 
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I  ?  3 4  5  6  7  8  9  0  1  2  3  4  5  6  "*  8  9  0  1  2  3  4  5  6  7  3  9  0 1 2  3  4  5  6  7  8  9  0  1  2  3  4  5  6  7  8  9  0 1 2  3  4  5  6  7  8  9  0 1. 2  3  4  5  6  7  8  9  0 1 2  3  4  5  6  7  6  9  u 

PRIORY07  1   2  28  48 

PRIORYOR  .  1  34  48 

PRIORY 09  4  8 

PR  I  OR YIO  49  4  8 

PRIORY11   1  48 

PRT0PY12   1  48 

PRIORY! 3   1  4  8 

PRI0RY14   1  48 

PRT0RY15   1  4  8 

PRIORY  16   1  5  0  48 

PRIORY17   1  48 

PRIORY18  52  5  3  5  4  bb    4b 

PR  I  CRY  19   1  5  6  4  \ 

pojopy?n   1  57  5r^  54  53  b2    48 

W  r  I C  A  P  0  1 

.'rLr  a  pn?  1  44  144  144  144  144  144  144  144  144  144  144  144 

WfrLCAP0313  13   13   13   13   13   13   13   13   13   13   13 

WELCAP040  5  05   05   05   05   05   05   05   05   05   u5   o5 

WELCAP0505  05   05   05   05   05   05   05   05   05   05   U5 

WFLCAP06072  0720^72  072  072  072  072  072  072  072  072  ul  Z 

,'"rLCAP07144  144  144  144  144  144  144  144  144  144  144  144 

'•"rLCAP0Rn86  0P6  086  036  086  ^36  036  086  086  086  086  066 

wcLCAP091  11111111111 

'■'rLCAP10 

W,rLCAP11144  14441*-  ♦■  1  (.4  144  144  144  144  144  144  144  144 

WELCAP125  5  :■  555555555 

WELCAP13 

WELCAP14065  065   )6 5  065  065  065  065  065  065  065  065  065 

1'.'ELCAP1504  04   0  \        Q'»   04   04   04   04   04   04   04   04 

WFLCAP16-072  072  072  0^2  072  072  072  072  072  072  072  072 

WFLCAP17 

,'."rLC./\P18065  065  06  5  065  065  065  06  5  065  065  065  065  065 

WELCAP1905  05   05   05   05   05   05   05   05   05   05   05   u5 

WELCAP20 

WFLCAP21 

VPLCAP??1?2  12?  I7?  122  12?  122  122  122  122  122  122  122 

WELCAP23108  108  103  1081108  103  108  106  106  lOd  i06  108 

W  E  L  C  A  P  24 

'/rLCAP25144  144  1.44  144  144  144  144  144  144  144  144  144 

WELCAP26094  094  094  094  094  0^4  094  094  094  094  094  094 

WELCAP27 

W  E  L  C  A  P  2  8  1  0  0 

WELCAP29 

WELCAP30 

WELCAP31144  144  144  144  144  144  144  144  144  144  144  144 

WE L C  A  P 3  2 

VELCAP3310S  108  103  108  108  103  108  106  108  106  108  108 

WELCAP34043  043  043  043  043  043  043  043  043  043  043  043 
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036 
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03  6 

WPLCAP36032 

032 

0-52 
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032 

032 

032 

032 

032 

032 

032 

032 

WELCAP37 

W  F  L  C  A  P  3  8  0  5 

05 

05 

05 

05 

05 

0  5 

05 

05 

05 

0  5 

05 

'•'FL^P3bl  1  5 

11  5 

':  15 

115 

115 

1  15 

115 

115 

1  15 

115 

115 

115  115  115 

W  E  L  C  A  P  4  0 

WELCAP4104 

04 

04 

04 

04 

04 

04 

04 

04 

04 

04 

04 

WELCAP421 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

WELCAP43 

WELCAP4405 

05 

05 

0.) 

05 

05 

05 

05 

05 

05 

05 

05 

WELCAP45 

WELCAP460-8  5 

08  5 

08  5 

08  5 

085 

085 

08  5 

085 

085 

Oo  3 

085 

083 

WFLCAP47029 

02^ 

)29 

Q-39 

02  9 

0  29 

029 

02  9 

02V 

029 

029 

029 

l*'ELCAP48014 

014 

J14 

0  \4 

014 

014 

014 

014 

0  14 

014 

014 

014  014 

WELCAP4905 

05 

05 

Oj 

05 

0  5 

05 

05 

05 

05 

0  5 

vb 

WE L CAP 50 

W  E  L  C  A  P  5  1 0  6  4 

064 

064 

064 

064 

064 

064 

064 

06 

064 

064 

u64 

'''rLCAP5  20  5 

05 

05 

0  5 

05 

0  5 

05 

05 

05 

05 

0  5 

0  5 

W  r  L  C  A  P  5  3 

!>'!"LCAP54 

WELCAP5  50  8  6 

086 

086 

08  6 

08o 

0o6 

086 

086 

086 

06  6 

08b 

utio 

W  E  L  C  A  P  5  6 

W  E  L  C  A  P  5  7  0  1  8 

013 

°18 

018 

018 

016 

018 

Old 

018 

0  i  6 

018 

uid  jIo  jib 

W  F  L  C  A  P  5  8 

'•'FLPCTC105 

06 

0  7 

02 

05 

11111111111103 

03 

03 

0  3   06        0  2   Oj 

WFLPCT02]  1  1 

I 

01 

01 

09 

04 

07 

111] 

1111          0  5   0  5 

WELPCT03 

WFLL   01 

WELL   0  2   1 

-> 

3  (- 

^   5 

11  ' 

12  13 

WELL   03   6 

7 

W ELL   04   8 

9  10 

WELL   05  14 

15 

WELL   06 

WELL   07  16 

17 

WELL   08  18 

19 

WFLL   09  20 

WFLL   10 

W  E  L  L   11  21 

. 

" 

WFLL  ■  12  22 

23  3?  34 

y 

WFLL   13  2  4 

25 

WFLL   14  2  6 

27 

WFLL   15  29 

WELL   16  2  8 

^0  31  3P 

i 

WELL   17  3  5 

36 

WELL   18 

W  f7  L  L   19  3  7 

WELL   20  38 

IMPMAX01 

03        1111 


51 


Appendix  3  -  Listing  of  the  Data  Cards 
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APPENDIX  IV 


SAMPLE  OUTPUT 


The  following  is  part  of  the  output  of  a  run  made  with  example  data  in  the  Ipswich  River  Basin. 
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Appendix  4  -  Portions  of  Model  Output 
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APPENDIX  V 


WELL  EFFECTS 


The  following  is  a  listing  of  the  Fortran  IV  source  deck  of  the  program  used  to  compute 
the  effects  of  well  pumpages  on  the  river;  a  list  of  the  data  cards  used:  and  part  of  the 
output  from  the  program. 
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Appendix  5  -  Well  Effect  Program,  Data  &  Output 


C 

C  PRCGRAM  FCR  CCMPUTING  WELL  EFFECT  ON  STREAMS 

C 

c 

C  SEE  U.S.  GEOLOGICAL  SURVEY  WATER  RESOURCE  INVEST  IGAT ICN 

C  SERIES  -  ECCK  4  CHAPTER  CI  BY  C.T.  JENKINS  FOR 

C  PHILOSOPHY  . 

C 

C  THIS  PRCGRAM  WAS  CCDEC  BY  I.C.  JAMES, II,  U. S .GEOLOGIC AL 

C  SURVEY,  WATER  RESOURCES  DIVISION,  3CST0N  MASSACHUSETTS. 

C 

c 

DIMENSION  C( 12)  ,  G0C12) 

C 

kPITE (6,167) 
C 

C      TRANS  =  TRANSMISSIVITY 
C      S  =  STORAGE  COEFFICIENT 
C      A  =  DISTANCE  FROM  STREAN  TC  WELL 
C 

5  REAC  (5,  1,END=2G0)  £ , TR AN S , S , NCWELL 

TRANS  -  TRANS/7. 5 

WRITE  (6,2)  NOWELL 
C 

CLP=C.O 

SN=C.C 

DC  10  1=1,  <5 

T  =  3.04*FL0AT(I  ) 

F  =  l.-ERF  ( A/SQRT(4.*T*TPANS/S) ) 
SN  =  cm  +  p 

10  CONTINUE 
T=30.4 

F  =  l.-ERF  (  A/SQRT(4.*T*TPANS/S)  ) 
SM  =  SM  +  F/2. 
Q(l  )  =  SM/10. 
CLM  =  CUM  +  C(  1) 
CC  20  J=2,12 
SN=0.C 

TL  =  20.4*FLCAT( J-l) 
SM  =  SM  +  F/2. 
DC  15  1  =  1,  c 

T  =  3.C4*FLOAT( I  )  +  TL 

F=ESF(A/SCRT(4.*(T-30.4)*TRANS/S)  )  -E  RF  (  A/  SQR  T  (  4  .  *T*TR  AN  S  /  S  )  ) 
SI*  =  SM  +  F 

15  CONTINUE 

T  =  3C.4*TL 

F=ERF(A/SCRT(4.*{T-3C.4)*TPANS/S)  )-ERF  (  A/  SORT  (4  .  *T*TRANS  /S  )  ) 

SM  =  SM  +  F/2. 

C(J)  =  SM/10. 

CLM  =  CUM  +  Q( J) 
20  CCNTINUE 

DC  21  J  =  1,12 
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Appendix  5  -  Well  Effect  Program,  Data  &  Output 


21    CC(J)    =    CUJ/CUM 
DC    25    J=13,12C 

T  =    30.4*FLCAT( J-l)     +     15.2 

F=EPF(A/SCRT(4.*(T-30.4)*TRANS/S)  )-ERF ( A/SQRT ( 4 . *T*TRANS /S  )  ) 

CUV  =  CUV  +  F 
25  CONTINUE 

WPITE(6,101)  A, S, TRANS 

OSUM  =  c. 

DC    30     1=1,12 
0(1)    -    Q(  I  )/CUV 
30    QSUy    =    QSLV    +    C(  I) 
*PITE(6,1C2)     Q 
lrfPITE(6f103)     CSUy 
WRITE<6,1C4)     CC 
GC    TO    5 
200     STCP 
C 
C 

167    FCPNAT11H1  ) 

1  FCP^AT  (3F10.CI5) 

2  FORMAT  (//////////1X, 'WELL  NUN9ER',I5) 

101  FCP^AT  (  1X,'CISTANCE  A=',F5.C,'  FT,    STORAGE  COEF=«,F5.3, 
$TRANSMSSIVlTY  =  »,F8.Ct//) 

102  FCPMATdX,  'MONTHLY  RESPONSE  ',  12F9  .  5  > 

103  FCPFAT(1X,'12  NONTH  CUMLLATIVE  RESPONSE= • , F9 . 5 ) 

104  FCRMATt  IX, 'RESCALED  RESPON SE •  ,  12F9 . 5  ,/// /  ) 
END 
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Appendix  5  -  Well  Effect  Program,  Data  A  Output 
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